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We present new sets of fragmentation functions in next-to-leading order QCD 
that are determined from e^e~ annihilation data of inclusive particle production. 
In addition to the 0{as) unpolarized cross section the longitudinal cross section 
is also used to extract the gluon fragmentation function from e~^e~ annihilation 
data. As the 0{a^) vanishes for longitudinal polarized photons (or Z bosons), 
the 0(0^) corrections are required to reduce the scale ambiguities. Recently, 
P.J. Rijken and W.L. van Neerven presented the longitudinal coefficient func- 
tions to next-to-leading order. We confirm part of their results in this thesis 
and complete the calculation by the results for the color class CpTji that must 
be included for a consistent comparison with LEPl data. The complete set of 
coefficient functions is then used together with novel data from ALEPH to deter- 
mine the fragmentation functions for charged hadrons. This set, and also sets for 
charged pious, kaons, and D* mesons as well as neutral kaons published previ- 
ously, can then be employed to test QCD in e^e~ annihilation, photoproduction, 
77 collisions, pp scattering and DIS. Finally, we suggest how the improved knowl- 
edge on the fragmentation in particular of the gluon could be used to determine 
the gluon and charm content of the photon. 




Abstract 



* Supported by Bundesministerium fiir Forschung und Technologie, Bonn, Germany, under Contract 
05 6 HH 93P (5), and by EEC Program Human Capital and Mobility through Network Physics at 
High Energy Colliders under Contract CHRX-CT93-0357 (DG12 COMA). 



1 Introduction 



At the heart of high-energy physics hes the desire to find the fundamental building 
blocks of matter and to understand how they make up our world. This is a very old 
quest that has made substantial progress only in the past 100 years with the advent of 
high-energy experiments, as the spatial resolution achieved in experiments is inversely 
proportional to the energy transfer. The picture that has emerged contains leptons (like 
the electron), quarks (that build up hadrons like the proton and gauge bosons (like 
the photon) that mediate the four fundamental forces gravitation, electromagnetism, 
weak, and strong interactions. Due to its weakness, gravitation can be disregarded in 
high-energy experiments. 

The first complete theory was formulated for electromagnetism. That theory. Quan- 
tum Electrodynamics (QED), succeeded in describing all experimental data on elec- 
tromagnetic interactions with a precision rivaled by no other theory to date. However, 
QED is not exactly solvable, it is solved in a perturbative expansion of its coupling 
constant oqed- As that is of the order of 10~^, neglecting terms of C(aQED) '^^ even 
C(aQE£,) introduces but a small uncertainty. Because of its tremendous success, the 
concept of a perturbatively solvable field theory was adopted to formulate theories for 
the other forces as well. In this process the weak interactions turned out to be closely 
related to the electromagnetic ones, leading to a unification of the two theories to the 
theory of electroweak interactions 0. 

The field theory for the strong interactions is Quantum Chromodynamics (QCD) 0]. 
Like QED it can be solved only perturbatively (the alternative lattice techniques are not 
applicable to the subject under study, at present, and are not discussed). In contrast to 
QED, however, the coupling constant as is much larger, and the fundamental particles 
of the theory, the quarks, are only free in the limit of asymptotically high energies 
This corresponds to small distances, i.e. the quarks are not free on macroscopic scales 
but are confined inside hadrons made up of two or three quarks in the case of mesons 
or baryons, respectively. 

Therefore, only part of the strong interactions can be calculated perturbatively, 
the confinement can not be computed from QCD (yet), and its properties must be 
measured. For a number of high-energy reactions, the hard and soft regions can be 
disentangled and treated separately. This is the content of the factorization theorem 
of QCD. The hard cross sections for parton scattering can then be calculated to fixed 
order in perturbation theory. Beyond leading order, virtual corrections to the process 
under study must be taken into account, too. These give rise to singularities in certain 
regions of the phase space, as do the real corrections. For inclusive processes, the 
infrared and collinear singularities cancel when the real and virtual contributions are 
added 0. The remaining ultraviolet singularities are removed by wave function and 
coupling constant renormalization. 
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For semi-inclusive processes this does not hold true anymore. In the case of single 
hadron production, e.g., the transition of the parton from a hard, quasi-free particle 
to confinement within the observed final-state hadron can not be calculated from 
perturbative QCD. It must be modeled by a phenomenological fragmentation function 
(FF). Because this function applies to a single parton, some (collinear) final-state 
singularities remain. These residual singularities are factored off into the bare, infinite 
FF, leading to finite partonic cross sections and finite fragmentation functions (an 
analogous procedure applies to initial-state singularities and parton density functions, 
PDF's). 

Inclusive particle production (IPP) is of particular interest for the investigation of 
QCD. It has been measured in a variety of processes and over a wide range of energies. 
Due to the well defined final-state, the measurements have small systematic errors and 
the observables are also unambiguously defined theoretically. This enables a test of 
QCD via two of its central properties: scaling violation and universality of factorization. 
As pointed out, the FF's required for the calculation of hadron production cross sections 
are not calculable, but must be fitted to data. Once this has been accomplished 
with the help of data of one specific process at one particular center-of-mass (CM) 
energy, the FF's so obtained can be applied to any process due to the universality of 
factorization. Their scale dependence is determined by the Altarelli-Parisi splitting 
functions. The determination of these FF's from data, and their application in a 
number of phenomenological studies are the subjects of this thesis. 

The choice process for the extraction of FF's from data is e^e~ annihilation. From 
the experimental viewpoint it is distinguished by an abundance of data that have been 
collected over the past fifteen years at energies from 3.6 to 91.2 GeV 0. In addition, 
the background is fully under control so that not only statistical but also systematical 
errors are very small. From the theoretical viewpoint, the greatest advantage of e^e~ 
annihilation is that no additional non-perturbative input is needed, avoiding possible 
bias in the determination of the FF's. Moreover, the Born approximation is of purely 
electroweak origin, i.e. it is precisely known, and any deviation from it is due to the 
strong interaction. Last but not least, the phase space integrations required for the 
calculation of single particle production are particularly easy for this process, owing 
to the simplicity of the initial state. This permits to perform the integration over 
phase space analytically (which is not possible, e.g., for pp scattering), speeding up the 
numerical evaluation of the cross section tremendously, which is crucial for fitting. 

To salvage the most of the predictive power of QCD it is desirable to use data from 
one specific process only in extracting the FF's. One disadvantage of e^e~ annihilation 
in this context is that gluons are not produced at LO, leading to a small contribution 
to the unpolarized cross section and consequently to little information on their frag- 
mentation. However, the gluons do contribute a major part of the longitudinal cross 
section. From a comparison of this measurement with theory one can determine the 
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gluon as well. Until recently, the longitudinal cross section has been computed to 
LO accuracy only - which is 0{as) - introducing large scale dependences. The NLO 
{0{al)) corrections are required to reduce this uncertainty. The calculation of the 
relevant coefficient functions has been presented by P.J. Rijken and W.L. van Neerven 
for the color classes Cp, CpiCp — l/SAi'c), CpNc, and NpTjiCp. In this work we 
con&m the results for the class NfTrCf and augment the results by the color class 
CfTr. The contribution from this class turns out to be numerically small. 

Subsequently, we use the complete set of NLO longitudinal coefficient functions 
together with the NLO transverse coefficient functions to determine a set of charged 
hadron FF's to NLO that feature a better constrained gluon. Although the transverse 
cross section has been calculated to NNLO 0, which is of order 0{a1) and hence 
formally of the same order as the longitudinal cross section to NLO, these higher 
orders are not included in the present analysis, as they would lead to inconsistency 
with the Altarelli-Parisi evolution kernels that are only known to NLO. For later use 
in consistent fits to deep inelastic scattering (DIS) and photoproduction data, the 
new fit is performed for distinct values of A. The range from A*^^^ = 150 to 400 MeV is 
covered in 50 MeV steps, our central choice is A*^"^) = 350 MeV. 

In addition, a number of studies undertaken in the past years in collaboration with 
B.A. Kniehl and G. Kramer will be reviewed in order to demonstrate the potential of 
inclusive particle production as a tool to test QCD in a variety of processes. 

Specifically, we ffited charged pion and charged kaon fragmentation functions to 
data from TPC |^ and to data from ALEPH [|10],|lT|. The novel data in |[T^ for the sum 
over charged hadrons distinguished between different quark flavors. Together with data 
on the gluon FF from OPAL [|T2], this enabled us to determine the FF's for individual 
flavors for the ffist time We also presented sets of neutral kaon [14| and D*"^ meson 
FF's Owing to the universality of factorization, these FF's are applicable to any 
process. We employed them in photoproduction, pp, and 77 scattering and found good 
agreement with data. The main emphasis was put on photoproduction, where very 
good data from HERA are available for study ||T6| , p!7| , p!8| , p!9| , pO| , ^ . In particular the 
rapidity spectra have the potential to improve our knowledge of the gluon and charm 
content of the photon [|15|,Q. In 77 scattering, NLO predictions for inclusive particle 
production at LEP2 and at the NLC were presented for the ffist time |23[] . 

The thesis is organized as follows. Section 2 contains a short review of the theoretical 
framework for using perturbative QCD and a discussion of the implications of the 
factorization theorem. In section 3 the process of e^e~ annihilation is discussed with 
emphasis on the longitudinal cross section. Sets of fragmentation functions for various 
mesons, extracted from e^e~ annihilation data, are presented in section 4, and with 
the help of the 0{a1) longitudinal coefficient functions a new set of particularly well 
constrained charged hadron FF's is obtained. Subsequently, these sets are applied 
to test QCD, via scaling violation and the property of universality, in section 5. A 
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new approach to constraining the (incalculable) soft regime of QCD in the case of the 
photonic gluon and charm densities is also presented. The results are summarized in 
section 6 together with an outlook on the potential of future analyses. 
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2 Factorization of Hard and Soft Regions 



The strong coupling constant that usually serves as the expansion parameter in pertur- 
bative QCD becomes small only for short distances. This corresponds to high momenta 
so that the bound states of colored objects can not be calculated and the problem of 
confinement remains unsolved. To put perturbation theory to any use in QCD, the 
short-distance and long-distance processes therefore have to be separated. Only for 
processes and experimental observables where this is feasible a comparison between 
perturbative QCD calculations and experimental data is possible. 

For this reason the factorization theorem of QCD that ensures the factorizability 
of hard and soft regions (corresponding to short and long distances, respectively) for 
a number of processes is crucial to any calculation in perturbative QCD. Moreover it 
is of particular interest for this analysis where the non-perturbative fragmentation of 
partons into hadrons is treated expressly by the parametrization of the ignorance about 
the confinement into phenomenological functions. 

In Fig. 0, the fragmentation process is represented by the large blob to the right. 
The double lines with the arrowheads represent hadrons in the initial or final state, 
the blob labeled H represents the calculable hard scattering and the small blobs in the 
initial state represent another non-perturbative ingredient required for some processes, 
namely parton density functions of hadron^. 




Figure 1: General structure of factorization of a high-energy cross section. 

This schematic diagram shows the general features of factorization in any factoriz- 
able process. The large blob to the right may, however, also be taken to represent jet 



fragmentation instead of inclusive particle production. The reader is referred to p4| , p5 
and references therein for details on the subject of jet production. 

The purpose of this section is to give a short account of the content of the fac- 
torization theorem and to discuss its consequences in some detail. The theorem will 
be employed in section 4 to fit FF's to data from e~^e~ annihilation and some of its 

^ Here as in all Feynman diagrams in this work, time flow is from the left towards the right. 
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implications will be used as rigorous tests of QCD when applying them in section 5. A 
more detailed intuitive picture of factorization can be found in ||26 |. 



2.1 Factorization Theorem 



The factorization theorem states that the soft and the hard parts of the cross section 
factorize (in fourier space) so that the cross section can be written as the convolution of 
a non-perturbative part with a perturbatively calculable hard scattering cross section. 
It holds in every renormalizable field theory for all leading twist two contributions, 
in particular, it has been proven in |^ that the inclusive cross section for any hard 
process factorizes. To leading order, both the partonic, hard cross section and the 
FF's are finite. In higher orders, however, singularities are generated through the 
real corrections. These are in part canceled by singularities in the virtual corrections 
that have to be added. Due to the Kinoshita-Lee-Nauenberg theorem, the cancellation 
works perfectly for inclusive quantities . For semi-inclusive observables, like inclusive 
particle production, residual singularities are factored off into the bare FF's. This 
is achieved with universal transition functions F. In the dimensional regularization 



scheme [28,29| the space-time dimension n is taken slightly off the physical 4 dimensions 
in intermediate steps of the calculation. This procedure isolates the singularities as 
poles in 4 — n = 2e. The transition functions have the general structure 
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Se = r(i - 

are artefacts of dimensional regularization. 



5) (471) 



(2.1) 



(2.2) 



(2.3) 



Mf defines the scale at which the diver- 
gences are subtracted. The unrenormalized strong coupling constant is denoted by TXs- 
Dependence of the functions on variables is suppressed to simplify notation; 1 is short- 
hand for the 5 function. The superscripts T, S stand for the timelike and the spacelike 
processes, respectively. The formulae are fully analogous and as e^e~ annihilation is 
timelike (Q^ > 0), only the timelike functions will be considered in the following. The 
superscript will be dropped for ease of notation. The da^b functions specify the fac- 
torization scheme, the MS-scheme adopted in this work is defined by da^b = for 
all a, b. 
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Up to order a^, the timelike F functions in the MS-scheme are given by: 



\ p(o) /2 4^ 



2 



e 



e2n ^ \Mf 



K^, = (2.5) 

= rfi, + rfi, (2.6) 
rG^. = -^^^^(^') (2.8) 



2 \ ^ 



The singularities in the bare cross section da are factored into the bare fragmentation 
functions D according to 







(2.10) 




= rJL^ (g) dalg + Tc^q (E) daL,G 


(2.11) 




= 2Tq^G ® daf^g + Tg^g ® daL,G 


(2.12) 



and analogous equations 

d';,s = K',®d''^s (2.13) 

= rj_., ® + 2NfT,^g ® (2.14) 

= Tg-^g ®Da + -^^G->, ® % (2.15) 

for the FF's. The definition of the combinations S and NS of FF's in the above 
equations will be given in the next section. The physical hadronic cross section remains 
unchanged by this operation, 



da'' = Y,dai(E = J2dai^D'^ . (2.16) 



_ —h 

^aai (X 

/ I 



The details of the summation are given in (|3.5|) . The shorthand (S> denotes convolution 

\fi^g]{x)= [ dxi I dx2f{xi)g{x2)5{x — X1X2) 
Jo Jo 

"1 dz 



(2,1.) 



The actual factorization for color class NfTrCf will be performed in subsection 3.3.3. 



One of the main results of the factorization theorem is that the functions parametriz- 
ing the non-perturbative part are process independent. This feature, known as uni- 
versality of the FF's and the PDF's, is instrumental to the predictive power of QCD 
analyses. It is one of the most directly testable properties of QCD and one of the rea- 
sons for the phenomenological significance of inclusive particle production and FF's. 



Albeit, there are some caveats to the predictive power of the factorization theorem 
as applied in this work. 

(i) As stated at the outset, the factorization theorem holds for leading twist two 
contributions, i.e. for terms of Oil/Q"^). Thus, the results of present day perturbative 
calculations are correct at least in the limit — > oo. For any finite energy, higher 
twist of the order l/Q'^, > 3 may be important, too. For deep inelastic scattering 
it has been proven that the subleading twist is twist four, hence those contributions 
are suppressed by a relative factor of and may safely be neglected. In e^e~ 

annihilation, the situation is less clear as the subleading terms may be of 0{1/Q^) 
for some observables. We demonstrated some time ago that present data on the 



unpolarized inclusive particle production do not necessitate the introduction of such 
terms as even the scaling violation, where they would show up most prominently, is in 
good agreement with data. This is demonstrated, for example, by Fig. Q in section 5.2. 
Very recently, M. Dasgupta and B.R. Webber |^ calculated the corrections with the 
help of renormalon techniques and found that the 0{1/Q^) is absent in unpolarized 
IPP cross sections also in e^e~ annihilation. In the same work they also maintain that 
the longitudinal and transverse contributions separately do have contributions. 
The uncalculable coefficient of those contributions, however, is small and omitting 
such terms altogether is consistent with the data on Ree,L Ol at the 1.5 a level. Power 
corrections will thus be neglected throughout this work. Some phenomenology of higher 
twist and power corrections in e^e~ annihilation can be found in P^p3| , pB[ . 

(ii) How much of the cross section is factored into the non-perturbative functions 
is unequivocal only when the calculation is performed to all order^. For a finite order 
calculation the scale Mf sets the point at which the infrared singularities are factored 
off, introducing an additional unphysical parameter in the results of the calculations. 
The main reason why one needs to calculate at least the subleading order for most 
practical purposes is to reduce the dependence of the result on this scale (and on 
related scales for the renormalization of the coupling, /i^, and for the factorization of 
the initial state singularities. Mi). In subleading order, called NLO in the following (i.e. 



^ Strictly speaking, this is not feasible even in principle, as the perturbative expansion in as is 
an asymptotic series, i.e. docs not converge. For the sake of the argument one may pretend it is a 
convergent series as only the first few orders will likely ever be calculable in practice, anyway. 
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NLO to a process that receives its leading contribution in 0{a^) is of order O^a^^^)), 
the sensitivity to the choice of scales is reduced compared to the LO. For actual fits to 
data it is quite cumbersome to vary the scales as one has to produce a separate fit for 
every choice of scales. For this reason the discussion of scales in this work is limited 
to the demonstration of the reduced scale dependence in NLO versus LO calculations. 
In the actual fits and applications the renormalization scales and fragmentation scales 
will generally all be identified with an appropriate hard scale of the process. In e'^e~ 
annihilation, there is only one hard scale, the virtuality of the intermediate vector 
boson, which somewhat reduces the ambiguity of choice of scales in this process. 

(iii) Last but not least, all particles, hadrons and partons alike, are treated as mass- 
less in this work. For the partons, the non-negligible masses of the c and b quarks are 
usually taken into account only via thresholds in the evolution of the FF's, accord- 
ing to the variable flavor number scheme (VFNS) 1^^. In the recent past, also some 
massive calculations of cross sections have been performed [^^]. In this thesis, the 
VFNS philosophy is adopted. Only for the heavy D*^ mesons this topic is considered 
in more detail and the factorization scheme is modified also to properly account for 
both the intermediate- and high-energy regions. As for the hadrons, there is no easy 
way to account for their masses. The masses are thus taken to vanish throughout. 
For momenta that are large compared to all the masses involved, this simplification is 
well justified. However, when the momentum of the final-state hadron is of the same 
order as its mass, threshold effects arise and the massless results become unreliable. 
This is a general feature of all massless calculations and is also observed, e.g. in jet 
production when approaching the boundary of phase space. In e^e~ annihilation the 
non-vanishing quark masses are important only for small virtualities of the photon. 
The finite hadron masses always become important at some point for production of a 
hadron with small momentum fraction x, notably for heavy hadrons like D*^ or pro- 
tons/antiprotons. On the other hand, the methods employed traditionally for fitting 
and evolving FF's become unreliable for very small x in any case. This is due to the 
large logarithms ofl/x that have to be resummed, eventually. The Altarelli-Parisi evo- 
lution thus becomes unreliable at some point and double-leading-logarithmic evolution 



equations are relevant instead |3^. In an intermediate region of moderately small x, 
the description of the fragmentation process can be improved by the modified lead- 
ing log approximation (MLLA) in connection with local parton hadron duality [^ ,^. 
This fundamental limitation of the traditional approach turns up as a singularity in 
differential NLO cross sections at x = 0. In this work, Xmin in the range 0.05 ... 0.20 
will be taken as lower limit for the validity of the massless approach. 
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2.2 Altarelli— Parisi Equations 



The fragmentation functions that have been introduced in the previous subsection are 
non-perturbative. They can not be calculated, at present, but must be taken from 
data at some scale fi. Their evolution to other scales, however, is then determined 
by perturbative QCD. In this subsection, we will present the equations governing the 
evolution of FF's. For brevity, they will be referred to as the Altarelli-Parisi (AP) 
equations, for a more complete bibliography of contributions see |^ . 

Let D^{x,iJ?), D^{x,fi'^), and Dq[x^ij?) be the FF's of Np quarks g, antiquarks 
g, and the gluon G, respectively, into some hadron h with momentum fraction x at 
the scale /i. The /x^ evolution of these FF's is conveniently formulated for the linear 
combinations 



Np 



1=1 



(2.18) 
(2.19) 
(2.20) 



as for these the gluon decouples from the non-singlet {+) and the asymmetric (-) com- 
binations, leaving only the singlet and the gluon fragmentation functions coupled; 



/i 



2 d 



d 
djj? 



I^f,),(x,/x^) = [P(+)(a,(/x2))®Df,),,(/x2)](x) 



Df_),,(x,/x2) = P(_)(a,(/x2))®Df_),,(/x2) (x) 



(2.21) 
(2.22) 



dij? 



+2Nf 



dfj? 



2Nf 



P< 



(a,(/.2))®D|(/i^ 



(2.23) 



+ 



Pr 



G-*G 



(a,(/i^)) D^^{^?)\ (x) 



(2.24) 



These equations can be derived from ( p.l3|) through (|2.15| ) by taking the derivative. 
The splitting functions that appear in above equations are defined as 



P(+) {x,as{n'^ 



q^q 



(x, as(/i^)) + {x, asif/)) 



(2.25) 
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(x, asifi^)) = P(+) (x, + ^NpPq^q (x, 

In the perturbative expansion of the sphtting functions, 



(2.26) 
(2.27) 



271 



(2.28) 



the first two orders have been calculated in . As the formulae are given in an implicit 
form only, in the x — 1 limit, we have collected the timelike functions in a ready to 
use form, i.e., with the coefficients of the delta functions and the plus distributions 
extracted explicitly, in |T^0. Here, we will list only the timelike LO splitting functions 
which enter the NLO cross sections via the scale dependences. 



Pflci-) 

Pi%i^) 
pSU^) 



2 VI — X 



1 



X 







1 + (1 -x) 



X 



2NfTr x2 + (1-x) 



^iVc - '^-NfTr ] 5(1 



X 



+2Nc 



— X/ + X 



2 + X - x^ 



(2.29) 
(2.30) 
(2.31) 
(2.32) 

(2.33) 



For flexibility in the parametrization at the starting scale, the evolution of the Da 
may be performed in x-space, solving the AP equations in their integral form 

D^^x, k) = D'^ix, 0) + r dk'^ Y: [Pa^i^Mk')) ® D^{k')] (x) (2.34) 

JO 27r ^ L J 

by iteration. The variable k = \n{fi'^/nl) replaces for ease of notation. The first 
term on the right-hand side of ( p.34| ) is the distribution at the starting scale /iq. 

Alternatively, the AP-equations can be solved analytically in moment space. Only 
the back-transformation has to be performed numerically in this approach. Details on 
the techniques of evolution in x- and n-space are given in appendix C. 



^The spacelike functions have been treated in a similar manner in 
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2.3 Momentum Sum Rule 



Any parton produced in the hard scattering process fragments with probabihty one 
into hadrons. In connection with local parton hadron duality this is the basis of 



the concept of jets where all hadrons originating from the same parton are grouped 
into a jet that is taken to bear a one-to-one correspondence to the primary parton, 
e.g., carries the same momentum. If one wants to study details of the jet one has 
to go beyond this and must investigate the probability of parton a to fragment into 
some specific hadron h. Due to momentum conservation, the momentum fractions 
/q^ dxxD^{x) of the individual hadron species must add up to one when summing over 
all hadrons h, 

r dxxD^^{x,Q'^) = l. (2.35) 



1 

h -0 

When fitting the related PDF's, e.g. of the proton, to data, the related sum rule 



I dxxFa/p{x,Q'^) 



is usually used as a constraint in the fit. Unfortunately, the sum rules ( p.35| ) can 
not be employed for fixing (some of) the parameters of the FF's for charged hadrons, 
because the charged hadrons, charged pions or other particles considered here form an 
incomplete subset of all hadrons and equation ( |2.35D is valid only for the sum over all 
hadrons. For the sum over a subset, the integral will give some value la for any parton 



a. These values could be extracted from the ALEPH |^5[ data — but to no avail as 
la 7^ h for two distinct types of partons a and h will, via mixing in the evolution, 
give rise to dependence of the la- As the FF's are parametrized at a lower scale 
yUo < no information on the relevant la values is available. 

In iTBIl, we used the fundamental momentum sum rule as an additional check of 



the FF's for charged pions and kaons. To account for the other hadrons, we had to 
make a number of assumptions. The protons and antiprotons were modeled to data at 
the Z pole by a function of the momentum fraction, as given in ( |4.4| ). Other charged 
hadrons were neglected. Of the neutral hadrons, only the pions and kaons were taken 
into account. Based on SU(3) symmetry, we could reasonably assume that they behave 
like their charged counterparts^. The inclusion of other neutral hadrons would have 
introduced unnecessary bias. When summing over these final states, values close to 
one were obtained. They exhibited a slight decrease with Q^, as shown in Table |l|. 
This is in agreement with the intuitive picture of the opening up of additional channels 
at higher energies. 



^For the pions, it has also been found experimentaUy | p6[ that the 7r° exhibits the same spectrum as 
the average of its charged partners. Only for very small x, where its momentum is not large compared 
to its mass any more, does their minute mass difference show up. 
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a 


Q [GeV] 


V2 


4 


10 


91 


200 


u 


0.81 (1.06) 


0.86 (1.01) 


0.86 (0.98) 


0.81 (0.89) 


0.79 (0.86) 


d 


0.93 (1.03) 


0.96 (0.99) 


0.94 (0.96) 


0.88 (0.88) 


0.86 (0.85) 


s 


1.03 (1.64) 


1.03 (1.50) 


1.01 (1.41) 


0.94 (1.25) 


0.92 (1.20) 


c 




0.90 (1.10) 


0.89 (1.05) 


0.84 (0.95) 


0.82 (0.92) 


b 






0.95 (1.12) 


0.90 (1.02) 


0.88 (0.99) 


G 


0.92 (1.07) 


0.96 (0.98) 


0.92 (0.91) 


0.81 (0.78) 


0.78 (0.75) 



Table 1: Left-hand side of eq. ( ^ISS] ) at NLO (LO) for Q = V2, 4, 10, 91, and 200 GeV, 
evaluated with the set |jl^. For technical reasons, the range of integration is restricted to 
X G (0.02,0.98). Hadrons are taken to be either pions, kaons or protons/antiprotons as 
specified above. 
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3 Particle Production in e+e Annihilation 



For the extraction of the fragmentation functions from data, the inclusive particle 
production (IPP) in e^e~ annihilation is the choice process. It has three big advantages: 

(i) It is particularly clean in the final state, i.e. no target or remnant jets contaminate 
the measurement. Hence, the background is fully under control, experimentally. 

(ii) No non-perturbative input is required for the comparison with QCD calculations 
(like the parton density functions in pp or ep scattering) , affording an accuracy in the 
determination of the FF's that is limited by the theoretical scale ambiguities (and the 
experimental errors), only. 

(iii) The process has been measured over a wide range of energies and, especially 
at the Z peak by LEPl, to a very high accuracy of the order of a few percent. 

For these reasons, this process has traditionally been employed to fit FF's (see 
section 4 for a review) and we will follow the tradition in this work. One disadvantage 
of e~^e~ annihilation in this context is, however that the gluon does not contribute to 
the cross section in LO, which leads to a poor correlation of the gluon fragmentation 
function and the data on unpolarized IPP. Some groups have therefore exploited data 
from other processes, e.g. pp hX, in addition to e'^e~ annihilation to constrain the 
gluon fragmentation function Dq p7| | . 

The philosophy of the analyses presented here is to salvage as much of QCD's predic- 
tive power as possible by using only the e~^e~ annihilation process in the determination 
of FF's. The property of universality can then be tested by applying the FF's to differ- 
ent processes as will be demonstrated in subsection 5.1. Following an idea by P. Nason 



and B.R. Webber |]3J], the longitudinal cross section o"/, is used to constrain the gluon 
FF instead in our new fit. For this part, the gluon and the quark both contribute with 
comparable strength to the IPP. (Tl vanishes in 0{a^) and receives contributions from 
the quarks as well as the gluons at 0{as). To reduce the scale dependence, one needs 
to include the NLO corrections also in this case, albeit of 0{a^^). 

In the following two subsections, the IPP process in e^e^ annihilation will be dis- 
cussed in general and the transversely polarized part will be briefly reviewed. In sub- 
section 3.3, we calculate part of the 0{al) corrections to the longitudinal process, 
amending and partly checking the results presented in The section closes with a 
short discussion of the new results. 
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3.1 The Process 



Consider the inclusive production of a hadron h m. e^e annihilation, 

e+(pa) + e-{p,) (7*, Z){q) ^ h{pu) + X . (3.1) 

Here, h denotes the measured outgoing hadron with four-momentum ph, the interme- 
diate vector boson is either a virtual photon or a Z with four-momentum q = Pa + Pb, 
and X stands for the sum of any other debris from the hard scattering as depicted in 
Fig. |. 




^ = -7^ , 0<x<l. (3.2) 



Figure 2: The IPP process in e+e annihilation. 

The virtual 7/^ is timelike, i.e., Q"^ = > 0. For large Q^, the transverse com- 
ponents of the momenta can be neglected. In the CM system, the outgoing hadron is 
characterized completely by its angle with the beam axis, 9, and the scaling variable 

2 pfe • q 

As the hadrons are considered to be massless, the variable defined above is the energy 
fraction of the outgoing hadron and does not coincide with the momentum fraction 
that is usually measured in experiments via the curvature of the particle trajectory in 
a magnetic field. For not too small x, the difference is negligible but for x not much 
larger than 2mh/\fQ^ threshold effects obviously become important and can not be 
treated properly within this framework. 

The unpolarized differential cross section of process ( p.l[ ) can be decomposed into 
several contributions with different tensor structures. For inclusive particle production. 



it is given by [HH 



- -(1 + cos^e)—^ + -sm^e—^ + -cose—^ , (3.3) 



dxdcosO 8 dx 4 dx 4: dx 

where the subscripts T and L denote the contributions due to transverse and longi- 
tudinal polarizations, respectively, and the asymmetric term, labeled A, is due to the 
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interference of the photon with the Z boson. Experimentally, the angular structure 
( p.3|) has first been confirmed by the TASSO collaboration Usually, experiments 
do not determine Q distributions for inclusive particle production measurements but 
rather integrate over the angles. This eliminates the asymmetric term in equation ( p.3|) . 
As the angular factors are normalized to their total integral, the single differential cross 
section is 



doc doc 



+ 



dal 
dx 



(3.4) 



The asymmetric contribution is therefore disregarded in this work and the reader is 
referred to for a discussion of a\. The transverse polarization state dominates the 
cross section, as demonstrated in Fig. The solid (dotted) lines show the ratios of 
the differential longitudinal to the unpolarized (T + L) cross section in NLO (LO). We 
will now turn to the calculation of the differential cross section. 



T3 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 08 0.9 1 



Figure 3: Relative importance of the longitudinal cross section. The solid (dotted) curves 
show the NLO (LO) results for daL/dax+L, evaluated at mz with our new set of FF's. 

The hadronic cross section is given by the convolution of the FF's with the partonic 
cross sections. 



1 da'i^ 
Nccfq dx 



c 



1 dz 



i=l 



Z 



1 da^n 



-,1V1J 



1 d^l 
Nc (To dz 



Np 



i=l 



1 

Nc<Jo dz 
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for polarizations k = T, L. The superscripts S and NS denote flavor singlet and 
flavor non-singlet combinations, as defined in section 2.2. Contributions from the class 
CfTr are labeled by F, suggestive of the fact that they vanish for photons due to 
Furry's theorem, see below. The cross section is normalized to the total cross section 
of e+e^ fJ'^^fJ'" for massless leptons. 



3Q2 



(3.6) 



where oqed is the electromagnetic coupling constant and its running is neglected. 
Nc = 3 counts the color states and takes care of the additional, unobserved degree of 
freedom. In fitting to data, the cross section will commonly be normalized to the total 
hadronic cross section, 

Nf 

^had = ReeNc E QqXQ^O, (3-7) 



i=l 



where the QCD corrections read to order al 



R, 



ee,L 



27r 




+iv.c.f-l^in(|)-nc(3) + lf; 



4 



+A'.r.C.(l„(|)+4C(3)-H^ 



Ree.T 
-Rpp 



Ree,L + Ree,T 



(3.8) 

(3.9) 
(3.10) 



Note that only the transverse polarization states contribute to 0{a^), whereas the 
0{as) corrections to the integrated cross section are solely due to the longitudinal 
polarization state. Ree,L has to be taken into account when studying the renormal- 
ization scale dependence of the longitudinal cross section. The terms of Ree,T are 
not needed in our fit. They can be found in |^ where eq. ( p.lO| ) has been used as 
an independent check on the NLO (respectively NNLO) longitudinal and transverse 
coefficient functions. 

The functions ^/) D^i^.M]) are the combinations defined in sub- 

section 2.2. In the applications considered in this work, only pairs of charge-conjugate 
particles are considered in the final state, i.e., h = + h~ (or + h for neu- 
tral hadrons). Charge-conjugation invariance of the strong interaction then ensures 
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D^,{x,Mj) = D^ {x,Mj) and D'^_^ -{x, Mj) vanishes. The summation index Np de- 
notes the number of (active) flavors in the VFNS. At NLO, Mj defines the scale where 
the divergence associated with the collinear radiation off parton a is to be subtracted, 
cf. subsection 2.1. 




0.5 r 



20 40 60 80 100 120 140 160 180 200 

Q [GeV] 



Figure 4: Dependence of the electroweak charges on the virtuality of the vector boson. 
The solid (dotted) lines show the ratios of Qup(Q^) (Qup(Q^)) over QdowniQ"^), respectively. 
At the Z pole, the relative weight of Qup(Q^) becomes maximal, whereas Qup(Q^) has its 
minimum there. 



The additional factors 



+ 



K)%(af)^ (^f)%(a?) 



Np 

1=1 



(3.11) 
(3.12) 



are due to the different couplings of the quarks to the vector bosons. The functional 
dependence on the CM energy is plotted in Fig. | for the ratios Qup(Q^)/Qdown(Q^) 
(sohd) and |Qup(Q^)|/Qdown(Q^) (dotted). For energies <C m| the propagator terms 
Pi{Q^) (for the interference of 7* and Z) and P2{Q^) (two Z propagators) given by^ 



1 



(m| - 



4 sin^ 9w cos^ 9\y 



[m. 



and 



(3.13) 



^ The expressions have been derived in the narrow width approximation. Deviations from improved 
formulae that take into account the energy dependence of the width [ pT] are sniaU compared to the 
strong interaction effects under consideration. 
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P2(g^) = 2 (3-14) 



are small and can be neglected. The usual couplings are recovered in that limit, 
QqAQ"^) ~^ ciiid the diagrams of class CfTr that give rise to (icr|^^ do not con- 
tribute to the cross section in accord with Furry's theorem [^. The complete set of 



vector and axial couplings of the photon and the Z boson, respectively, is given by 

v1 = -l a7 = 

"^up ~ +1 fermion i 

7 1 
^down 3 



f f = Ta^e - 2v2 sin^ 6w af = Ta^e 

= T^^q - 2vl sin^ 9w af = T^,q 



where q stands for either an up or a down type quark and where 6w is the Weinberg 
angle. The third components of the strong isospin are 





1 


^3,6 — 


~2 ' 


^3,up = 


1 


-^3, down 


1 

^2 ■ 



The bare, unrenormalized partonic cross sections da^/dx are calculated via 
d4l = ^We-'e- (7*, Z) a,...a,)\'FS^^^ 

^ ^ -L^.iJ^'^PSW , (3.15) 



2Q2 4Q4 '"^ 

where PS*-'^ denotes the phase space of an /-particle final state ai...a;. The flux factor is 
1/(2Q2)^ and 1/(4Q^) accounts for the average over initial spins and for the two photon 
(or Z) propagators. The ^-dependence of the lepton tensor is factored out according 
to ( |3.3| ) so that instead of the full lepton tensor, which gives the sum of (Jt + ctl, just 

V = 4e^?E^^4'^ (3-16) 
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must be contracted with the hadron tensor in eq. ( |3.15| ). The partonic cross section 
can then be written in the form 



a 



(2) 



L A 



(3.17) 



The dimensional regularization scheme |2^j2^ is used to keep track of the singularities 
in intermediate steps of the calculation inn = 4: — 2e dimensions. The modified entities 
in (|3.17|) are given byQ, 

'47r/i2y T{2-e) 



a 



(2) 



O"0 



g2 ) r(2 - 2e) ' 
Nc(l-e) 



-H 



fJ,U 



with 



2e 



PS 



(0 



A, 



.^2.(Z-2)pg(0 



(3.18) 
(3.19) 
(3.20) 



The mass parameter /x serves to keep the unrenormalized dimensionless also in n 7^ 4 
dimensions. The factor 
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ri-e) 



T{2-2e) \Q\ 



(3.21) 



will appear in the expressions for the Z-particle phase spaces and will thus cancel in 



PS^ For k = T the sum in ( p.l7|) runs over the two transverse polarizations A = 1, 2 
and for k = L just A = contributes. In the following, we define 



L A 



(3.22) 



where the couplings have been factored out of the hadronic matrix elements in ( |3.19|) . 
The factor 1/Nc will be absorbed in the color factors, and the additional factor 1/(1 — 
e) arises from the normalization to the Born-level quark-ant iquark cross section in 
n dimensions, a^'^\ following the notation of The factor cr*^^^ is universal, i.e., 

independent of the number of particles in the final state. Since the final result will be 
obtained hj e 0, e might be as well set to zero right away in this common factor 
and cr^^'' then cancels with the normalization factor in ( |3.5| ). The same holds for the 
factor 1/(1 — e), leading to 



1 rfa(') 



k,a 



dx 



(47ra,)'-2|M|2ApS^'^ 
dx 



(3.23) 



^ As usual, the purely leptonic vertex is calculated in 4 dimensions from the start so that the 
electromagnetic couplings give a factor of — (47raQED)(47rQ;QEDM^^)- 
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The renormalization of Ug, and the factorization of the poles in ( p.23| ) into the (singular) 
bare FF's will be performed in subsection 3.3.3. This will introduce the renormalization 
scale fifi and the factorization scale Mf. 



3.2 Transverse Cross Section to 0{as) 

The polarization tensor for the transversely polarized photon reads 



A=l,2 



(A)^(A) _ _ _ .(0) (0) 



(3.24) 



in the Feynman gauge. The longitudinal polarization tensor that is subtracted in 
( |3.24| ),is given in the next subsection. 

The right hand side of (|3.23|) is expressed in the form of the coefficient functions 
Ci™) of parton a in m-th order. The leading order does not receive virtual corrections 
and no singularities need to be factored off so that the coefficient functions are in 
that case just the real matrix elements. For higher orders, things are slightly more 
involved and will be discussed for the longitudinal case. Here, just the final result in 
is given. 



the MS-scheme 30 
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dz 


1 




Ncctq 


dz 


1 




Ncao 


dz 


1 


daT,G 


Ncao 


dz 



c 



NS,{0) 



2n 



p(0) ( 



z) In 



0, 
0, 
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2P, 



Ml 



, (n(l) 



(3.25) 

(3.26) 
(3.27) 
(3.28) 



The strong coupling is evaluated at the renormalization scale fiR- For the choice /i^ = 
Mj = Q^, the logarithms in (|3.25| ), (|3.28|) vanish and the coefficient functions alone 
determine the cross sections. The above expressions are given to NLO which is 0{as) 
in the transverse case. The dependence on the factorization scale is reduced at NLO. 



Furthermore, the K factor for the NNLO, 



K 



NNLO 



Observable 



NNLO 



Observable^^^ 



(3.29) 
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is not expected to differ significantly from 1. For example, the transverse part of the 



total cross section to NLO, evaluated with as{mz) = 0.12 is from equations ( |3.8|) 

through ( p:TU|) 

^NLO I 

Reef = = T-1I7 = 0-963 . (3.30) 

This compares favorably to [^] 

i?^^P^ = 0.943 ± 0.005 . (3.31) 



The coefficient functions have been calculated to 0{a1) |^ but the NLO expressions 
are sufficient for our purposes and have the advantage that the splitting functions can 
be implemented consistently at this order. 



3.3 Longitudinal Cross Section to 0{a 



I) 

As can be seen from ( |3.3| ), the longitudinal polarization state of the intermediate vector 
boson is selected by the angular distribution of the hadron produced in the inclusive 
process. In the framework considered in this thesis, the three-momentum of the hadron 
is aligned with that of its parent parton. Hence, the momentum of the parton produced 
in the hard scattering is linked to the gauge boson polarization. The fragmenting parton 
is labeled 3. The polarization vector reads 

ef = P^-^ (3.32) 
and due to current conservation the polarization tensor is then effectively given by 

,(0) (0) _ Q2 PSf^PSu _ 4 

At (!?(a°), only the transversely polarized photon and Z boson contribute to the 
cross section. To 0{as), the longitudinal cross sections of the subprocesses are given 
by 

' *^-"^<f)(.), (3.34) 



Ncao dz 271 ^''^ 
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1 

Nc^o dz 

Nc (Jo dz 
1 da^G 



Ncctq dz 2tt 
Inserting this together with the coefficient functions from appendix D in 



(3.35) 

(3.36) 

(3.37) 
leads to 



1 d(^L nr ^ 



Nqcto dx 



+ 
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Np 



i=l 




+ 0{al) . (3.38) 



Since eq. ( |3.38| ) is only given to LO in a^, the value here is not expected to be the 
same as the one in the NLO relations, i.e., when the NLO expressions are inserted in 
(|3.5|) . ALEPH and OPAL have recently presented data on dai/dx. In Fig. H, 
eq. ( p.38| ) is evaluated with the NLO set of [jl3[ and two-loop and compared to the 



ALEPH data |^5[. The result obtained (full curve) falls short of the data by a factor of 
two. At this point, one should keep in mind that eq. ( |3.38D is a LO prediction, and one 
should be prepared to allow for a K factor, which is typically larger than one. This 
K factor may be simulated by changing the renormalization and factorization scales. 
When choosing low scales of iir = Mf = 20 GeV, satisfactory agreement with the data 
is found. 

In fact, already the total QCD corrections, Ree,L = 1 ~ Ree,T, fall well short of the 
experimental data in LO, 



Re^L = 0.037 



, whereas 



0.057 ±0.005 [33] 



(3.39) 



The inclusion of the subleading corrections to the longitudinal cross section re- 
duces the scale dependence, as has already been argued. Moreover, only at NLO do 
the coefficient functions acquire a scheme dependence. Consistently using the NLO 
Altarelli-Parisi kernels with the NLO longitudinal cross section leads to a partial can- 
cellation of this scheme dependence and thus to a better perturbative stability. Our 
policy therefore is to calculate the longitudinal cross section to NLO when comparing 
with longitudinally polarized particle production data and to calculate the unpolarized 
cross section, which is the sum of the longitudinal and the transverse parts, also to 
NLO, when comparing to unpolarized data. In the latter case, the NLO contribution 
to the transverse part is of 0{as) so that the longitudinal part is consistently added 
at the same order of a^, although this is the leading contribution. 
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Figure 5: Scale dependence of the longitudinal cross section in e^e~ annihilation to LO. 
While the LO result with Mj = niz (solid) falls short of the ALEPH data |45] by a significant 
amount, a very low scale of Mf = 20 GeV (dotted) leads to satisfactory agreement. The plot 
is taken from [^], where the calculation with the leading order coefficient functions was 
misleadingly labeled NLO. 



The functions and the formula for the cross sections to 0{a1) have been 
presented in for the process e~^e~ — » 7* — hX. In the next section, we will 
complete the calculation with the class CfTr that contributes to e+e" ^ Z ^ hX in 
the case of an odd number of active flavors. For a cross check we will also repeat the 
calculation of class NpTjiCp- 

The two classes treated here belong to the process /Z — > qqqq that does not 
receive virtual corrections to order 0{a1) - this order is the LO for this final stateQ. 
Hence, only the real corrections must be calculated in the next subsections. 



3.3.1 Matrix Elements for the Real Corrections 

In this section all matrix elements are given for the process 7*/^' — >■ qqqq to 0{a1) in 
the dimensional regularization schemeQ. Those for the process 7*/^ qqGG can be 

In addition, the two-loop virtual corrections do not contribute to (Tl at all, as one of the vector 
bosons couples directly to the observed quark/antiquark in the relevant cut-diagrams, leading to 
vanishing matrix elements by current conservation. 

® The calculation of the gamma traces has been performed with the help of REDUCE . For color 
class CpTfi, the traces involve 75-matrices that have been implemented in [p3| for four dimensions 
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found in 




1 



2 



Figure 6: Born-amplitude for the inclusive production of a parton / in j* /Z IX. 

The Born approximation to the e^e~ annihilation process is given in Fig. As has 
been pointed out, it is of order Moreover, it does not contribute to the longitudinal 
cross section. To 0{as), the amplitudes in Fig. ^ contribute to the cross section, with 
color factor Cp- As this is the leading order, there are no virtual corrections. 




Figure 7: Amplitudes contributing to 0{as)- 

For the factorization later on one also needs the 0{e) corrections to the matrix 
elements of the three particle final states that are given in the literature to 0{e^), only. 
In n dimensions, the matrix elements are - apart from the color factor Cp ^ given by 



Q22 = (3.40) 
Q21 = (3.41) 

Qll = 16(l-e)^ (3.42) 



for the fragmentation of parton 1. As this is a quark (or antiquark), the elements are named 
Qkl for the interference of amplitudes k and / in Fig. |^. The momentum of the fragmenting 
parton defines the azimuthal angle so that the results for the matrix elements differ if some 
other parton is selected. When parton 3, a gluon, fragments, the matrix elements are named 
G22 and so on and read 

G22 = (3.43) 

only. This is sufficient for the present analysis, as no poles appear in class CfTr. One may obtain the 
n-dimensional results with the use of the Breitenlohner-Maison scheme [28,|4j for 75 that has been 
implemented in p5| , ^ . 



26 



G21 = 32^ 
Gil = 



(3.44) 
(3.45) 



Xi = 2pi-q/Q'^ is just the fraction of the total energy that is carried by parton i. The 
Uij are normahzed dot products of the momenta of particles i and j, c.f. ( p.88|) below. 

To order «g, many amplitudes contribute to the real corrections to the longitudinal 
cross section. Following the conventions of |0 , they are shown in Fig. ^ for the qqGG 
and in Fig. &) for the qqqq final-states, respectively. 



B 











a) 



b) 



Figure 8: Amplitudes contributing to ©(a^), a) for the qqGG and b) for the qqqq final state. 
Parton 1 (and 3 for b)) is a quark and 2 (and 4) an antiquark. 

When the amplitudes of the process |^ are added and squared, the interference (and 
diagonal) terms come with three distinct color factors. Those of process Hb also come 
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with three different factors. Following [Q, they will be grouped together according to 
their color factors, giving 6 contributions to the longitudinal cross section, labeled A 
through F, as hsted in Tables ^ and ^. 



perm, of 
first row: 


class A 


class B 
Cf{Cf - iiVc) 


class C 
CfNc 


(34) 
(12) 
(12)(34) 


All A32 A21 A22 
A44 A65 A54 A55 
A64 A66 
A31 A33 


A42 A52 A53 A41 
A51 A62 
A61 

A43 A63 


A71 A72 A82 A77 A87 
A74 A75 A85 
A84 A86 A76 A88 
A81 A83 A73 



Table 2: The color classes that contribute to the process e~^e qqGG. 



perm, of 
first row: 


class D 
NfTrCf 


class E 
Cf{Cf - iiVc) 


class F 
CfTr 


(24) 
(13) 
(24)(13) 


B77 B88 B87 
B55 B66 B65 
B33 B44 B43 
Bll B22 B21 


B83 B76 B73 B86 B84 B75 
B61 B85 B51 B62 
B74 B32 B42 B31 
B52 B41 


B81 B82 B53 
B63 B64 B71 
B54 
B72 



Table 3: The color classes that contribute to the process e'^e — > qqqq. 



The Tables are organized in columns of matrix elements that can be transformed 
into each other by interchanging a pair of partons in the final state, (12) for example 
represents the exchange of partons 1 and 2. Only part of the permutations in the 
tables are actually allowed for the longitudinal matrix elements, depending on which 
parton is singled out in the final state. For the fragmentation of parton 3, only the 
permutations (12), respectively (24) are permitted, whereas for the fragmentation of 
parton 1 the permutations (34) and (13), respectively, are applicable. 

The color factors are obtained from the color traces 



a 


= NcCf 


(3.46) 


tr {TaTh} 




(3.47) 


ifabctr {T^TfoTj 


= ^N'^Cf 


(3.48) 


tr {TaTbTaTb} 


= NcCf (Cf - ^Nc^ 


(3.49) 



28 



tr {TanT.Ta} = NcCl 



(3.50) 



In the case of the color gauge group SU(Ai'c), the Casimir operator of the fundamental 
representation has the eigenvalue Cp = [N^ — 1) / (A^c) and the trace normalization 

is Tr = 1/2 ig. 

For the process e~^e^ —>■ qqGG the complete set of NLO coefficient functions is given 
in 1^. We will hence only be concerned with the process e~^e~ QQQQ in the following. 
This process does not receive virtual corrections. We will recalculate the color class D 
as an independent check of the results presented in ^ and also consecutively calculate 
class F that is needed for the comparison to LEPl data. 





B87 



Figure 9: Cut-diagrams for e'^e~ annihilation at ©(q^), color class NpT^Cp- The cuts are 
not indicated - they run vertically through the center of each diagram so that the final state 
consists of four quarks. The fragmenting quark (3) is indicated by an X. 

Fig. § shows the cut-diagrams of class D. Any other matrix element can be obtained 
from those topologies via crossing, as specified in eq. ( |3.51| )ff. The fragmenting quark, 
marked by an X, couples only indirectly to the vector boson in the cut-diagrams 588 
and 587. The intermediate gluon is a fiavor singlet so that those matrix elements give 
rise to pure singlet contributions. These will be denoted by PS. Together with the 
non-singlet NS they make up the singlet, S. The matrix elements for class NpT^Cp 
are - apart from the common color factor - 



32 (2yi3 + 2^231/34 + Vli - £^34)^132/24 
^3^2342/34 



^gg _ ^^^23 ^ ^yiiy-j^: -r c<34 ^y34jyL:iy24: (351) 
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B87 = 


32 (ey|4 + 2yi3y23 + 2/132/34 + y23y34) (2/122/34 - yi32/24 - 2/142/23) 


(3.52) 




2^12/134 2/2342/14 


BU = 


32 ivh + yh 


- eyfr^ - 2eyi3y23 - ^2/23) 2/34 


(3.53) 




9 

^32/1232/12 


B33 = 


cr 24(511) 




(3.54) 


B77 = 


crl2(588) 




(3.55) 


566 = 


cr24(588) 




(3.56) 


565 = 


cr24(587) 




(3.57) 


555 = 


n(i42)(588) 




(3.58) 


544 = 







(3.59) 


543 = 







(3.60) 


522 = 







(3.61) 


521 = 


. 




(3.62) 



The operations cr 12 (M) etc. denote the interchange of two partons in the matrix 
element M. The 11(142) in eqn. ( 3.58| ) represents the cychc permutation of the indices 



1,4,2 in the manner 1 4, 4 ^ 2, and 2 — 1. The matrix elements for the class 
Cf{Cf — l/2A^c) are given by 



64 (2^23 + 22/232/34 + 2/34 - £ 2/34)2/132/24 

2/232/2342/34 

584 = (3.64) 



- Vyi 12— — 2— (3.63j 



^^'^ = ^ (e2/l22/l32/34 + ^^2/122/34 + ^yi3y2Ay34: " ^2/132/142/232/34 + £2/132/242/34 



-£2/142/232/34 - 2/122/132/232/34 " 2/122/132/34 + 2/122/232/34 " 2/132/232/24 
-2/132/242/34 + 2/132/142/23 + 2/132/142/232/34 " 3^132/232/242/34 

-22/132/242/34 - 2/142/232/34 ) ^ (3.65) 

^ 3^32/1342/142/2342/34 



„„„ 32 (£2/34 22/132/23 2/132/34 + 2/232/34) (2/122/34 - 2/132/24 - 2/142/23) 
576 = -2 2 (3-66) 

Q 3^32/1342/232/2342/34 
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575 = -7^ ( £2/122/34 - £2/132/242/34 + £2/142/232/34 + 2/122/132/34 - 42/132/14 



-2/132/24 - 32/132/142/23 - 22/132/142/34 - 22/142/232/34 ^ (3.67) 

^ a;§2/1232/1342/34 

574 = (3.68) 
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'^'^3 = (£2/122/132/34 -^£2/122/34 - £2/132/242/34 + £2/132/142/232/34 - £2/132/242/34 



+£2/142/232/34 + 2/122/132/34 + 2/122/132/34 " 2/132/24 + 2/132/142/23 

-2/132/242/34 - 32/132/142/232/34 " 22/142/232/34 ) (3-69) 

^ a;^2/i342/i42/34 
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32 / 2 

B31 = --^ (^2eyi3yi4 + 2eyi3yi4y23 + 2eyi3yi4y34 + 2e2/i42/23y34 



1-^12^13^34 - yuVsA - 2y?3yi4 - yhy24 

y23 



-yi3yi4y23 + ^13^24^34 - ^14^23^34) ^ (3.70) 

a;3yi22/i23yi34yi4 

562 = (3.71) 

561 = crl2(S83) (3.72) 

552 = (3.73) 

551 = crl2(573) (3.74) 

542 = (3.75) 

541 = (3.76) 

532 = (3.77) 

585 = crl2(576) (3.78) 

and those of class CfTr are 



DQ1 32 (7/14^23 - yi2y34 - yi3y24)(yi3y23 + 2/13^34 + 2/23) 

581 = —2 2 ^<-^(ej (3-79) 

Q 2^3^12^123^2342/34 

582 = (3.80) 

553 = n(i42)(581) (3.81) 

563 = cr24(581) (3.82) 

564 = (3.83) 

571 = crl2(581) (3.84) 

554 = (3.85) 

572 = , (3.86) 

where Bij stands for the interference term of amphtude i with amphtude j of Fig. ^jb. 
To account for the Bij, i < j that are not hsted, the off-diagonal matrix elements are 
weighted by a factor of two. Note that the matrix elements listed above are for the 
fragmentation of parton 3. As there are two quarks/antiquarks in the final state, the 
cross section receives an additional factor of 2. 



In Fig. |10| the cut-diagrams of the color class F are displayed. They differ from 
the cut-diagrams of any of the other classes (only partially given here, see for a 
complete set) in the number of gauge boson couplings per loop. While any fermion loop 
in the other classes couples to an even number of bosons, the loops 1 and 2 in class F 
each couple to three bosons, namely two gluons and one electroweak gauge boson Vi 
and V2, respectively. For the case of Vi = V2 = 7, the sum over the matrix elements 
vanishes due to Furry's theorem This cancellation comes about because there is 
essentially only one matrix element, 581, and the others can be obtained via crossing. 
However, permuting the momenta 1 and 2 in B71 gives a relative minus sign to 581 
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B81 



B71 



Figure 10: Cut-diagrams for e+e annihilation at ©(a^), color class CfTr. 

so that the two cancel, and analogously for the pair B63 and B53. When considering 
just one of the two vector bosons to be a Z-boson the same cancellation is at work 
and only the diagrams with Vi = = Z survive. Due to the axial part of the Zqq 
coupling (the vector part vanishes in the same manner as outlined above), the matrix 
elements are different from those with = 7. The main change is that they acquire 
an additional minus sign when crossed which spoils the cancellation. Apart from the 
minus sign, though, the symmetry is still intact and thus only 581 needs to be kept 
with an extra factor of four. As usual, parton 3 is taken to fragment. The electroweak 
coupling is then given by 

where the factor in square brackets accounts for the coupling on the leptonic side. As qi 
is unobserved one can right away sum over all g^, i = l...Np. Because of Gq^^ = —dq^^^^, 
the up- and down-type unobserved loops cancel. This would lead to a total cancellation 
if the quarks were truly massless. Albeit that is not the case, and at LEPl the top 
quark is not active in the VFNS framework, but the bottom must be taken into account. 
So the one graph that survives these numerous cancellations is 4581 with two Z's and 
one b loop. 

The various interference terms in each color class are finally grouped together. The 
matrix element for any specific color class then reads 

|M|2 = NcF^oior E BtJ , (3.87) 

i>j 

e.g., \M\l = CfTr{4:B81) for the class F. 

The Xi, Uij and yijk that appear in the matrix elements are simple functions of dot- 
products of the four-momenta of the intermediate vector boson and the four final-state 
particles, q and pj, respectively (i = 1, ..,4): 

Xi = , Xi e (0, 1] and = 2 (3.88) 

V i 
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Vrj = = g2 ' ^ 1°' 1] Z^^^i = ^ (3.89) 

...lol 

These definitions apply to any final state of / massless partons. In addition, for any 

specific / there are a number of useful relations between these variables that read for 
1 = 4: 

Uij = 1 - Xk - Xm + Vkm (3.91) 

Vij + Vik + yim = 'i-- Ujkm (3.92) 

Vijk — ^ (3.93) 

Vijk + Vijm + yjkm = 2 , (3.94) 

where k, m represent any permutation of 1, 2, 3, 4. 

For a cross check, also the unpolarized results (i.e. where the hadron tensor is 
contracted with —g^u) have been calculated to all orders in e and perfect agreement 
has been found with independent calculations in n dimensions . The n dimensional 
results for the classes A through E have already been presented in [p^. 

In the next two subsections, we will perform the phase space integration over the 
matrix elements of class D as a cross check and we will complete the set of coefficient 

F (2) 

functions by the calculation of the C^'^ . 



3.3.2 Phase Space Integration 



In this subsection, the integration over the degrees of freedom of the unobserved final- 
state particles will be performed. To this end, the phase space factors for /c-particle 
final states are given for k = 2. ..4 in the center-of-mass system of partons 1 and 2: 



(3.95) 



PS(^) = An 



647r2 



rfi 



An 



j dyi2dy23yiiy2iil - I/12 - I/23) ^©(1 - I/12 - I/23) 



(3.96) 



PS(^) = Ao 



16 647r4 



Stat 



1 I An 
T{l-2e) [q^^ 



2e 



[1 - 2e) 
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J dyi24dyi23dyi2iyi23yi24 - yi2) ^(yi2 + 1 - yi24 - yi23) ^2/l2 

0(2/l2)0(yi232/l24 - 2/12)0(1/12 + 1 - 2/124 " Z/m) 
"1 



, (3.97) 

where Aq is defined in (p.21| ). The expressions for PS*^^-* and PS*-^-* have also been derived 
independently and perfect agreement with has been achieved in n dimensions. The 



four-particle final state is determined by three invariants and the two angles 6 and 
(f). The angular integration over 6 is written in terms of f = 1/2(1 — cos^). The 
normalization factors for the angular integrations in PS*^^^ are given in appendix B and 
the statistical factor iVstat is due to the quantum statistics of the final-state particles. 
For the qqGG final state, A'stat = 1/2 and for the qqqq final state, Nst^t = 1/4. The 
two- and three-particle phase spaces are relatively simple and only the four-particle 
phase space that is the central ingredient in the calculation of the NLO coefficient 
functions requires some further attention. 

The coordinate system can be chosen so that, of the irreducible set of five invari- 
ants X3, 1/123, 2/12, 2/234, and j/23, only j/23 depends on cf) and only the last two depend 
on V. The angular integrations over the matrix elements 588 and Bll can then eas- 
ily be performed, the required integrals are given in appendix B. From (|3.94|) and 
( |3.93| ), 1/134 = 1 + X3 — ?/i23 — 2/234 and the denominator of 587 decomposes into inverse 
monomials of 7/234, 



1 1/11 



2/1342/234 1 + 2:3 - 2/123 V 2/134 2/234, 



(3.98) 



In the same manner, also the denominator of B81 gets decomposed. Moreover, in the 
case of B87, the matrix element is symmetric under (1^2), and so, by construction, 
is PS*-^'' so that the two terms arising from ( |3.98| ) can be treated collectively as 



64 (e?/34 + 27/137/23 + 2/132/34 + 2/232/34) ("2/122/34 + 2/132/24 + 2/142/23) ^^^ 
Q + X3 - 2/123)2/2342/34 

The two remaining energy integrations over 7/123 and 7/12 are restricted by the kine- 
matics (see ( |3.91| ) ff.), the boundaries are specified by the Heaviside functions in (|3.97| ). 



The resulting triangle is decomposed into two distinct regions, I and /J, as sketched 



in Fig. |ll| These regions are treated separately and the results are then added. The 



respective integration measures can be written as 

I,[/(x3,r,^)] ^ xl-'^{l-X3f-'' f'drr'-'%l 

Jo 



f 

Jo 



' dzz-%1 - zy ( 1 + (^ + ^3)r \ ^^^^^ ^ 
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yi23 



Figure 11: Schematic drawing of the four particle phase space. The shaded areas are the 

kinematically allowed regions. In the hatched region (I) yi2 is the singular variable, whereas 
for the cross-hatched region (/I) the combination j/123 — X3 — y\2 gives rise to poles in e. 



where the remaining invariants in terms of the new variables are 

yi2 = (1 - Xri)x2,rz and 
yi23 = x^r 

for region I. In an analogous way 

Jo 

[I - r)xs 



(3.101) 
(3.102) 



dzz-%l-zy' 1 + 
\ r 



z\ f{xs,r,z), (3.103) 



with 



1 H z j and 

1/123 = 1 - (1 - a;3)(l - r) . 



(3.104) 
(3.105) 
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From Fig. ^ it is obvious that in region / only yi2 and in region // only {yus — 
^3 ~ Uu) give rise to singularities in (|3.971 ). From the definition of the new variables 
in ( p.l01|) , ( p.l02| ) and ( |3.104]) , ( |3.105|) , respectively, one sees that in both regions the 
possible singularities then turn up in the z integration as 



dz z 



{1 + eg{z,X3,r)) 



(3.106) 



where g{z,X3,r) for z 0. The solution of the energy integrals both on and 
off singularities is briefiy sketched in appendix B. The last integration, of course, cor- 
responds to that over the momentum fraction x of the outgoing parton and is not 
performed. Information on the chosen coordinate system is obliterated by the phase 
space integration and the subscript of x will be dropped in the following. 

The integration over the matrix elements of class D yields a nonsinglet contribution 
(from Bll) and a pure singlet contribution. Of course, this is not the only class to 
contribute to the nonsinglet coefficient functions. However, the factorization works 
separately for each color factor so that it is possible to project out the terms propor- 
tional to NpTrCf and check the calculation for this part, independently. Together 
with the integration over the 0{as) matrix elements to order e, we obtain 



Nqctq dx 



yM^] 




NS,{2). 
L,q I 



X] 



(3.107) 



Nc<Jo dx 



1 daL,G 
Ncctq dx 



27r; ^ 2£ 



41. ® ^L,G 



.(1) 



fx) + C 



L,q I 



1 



,2 \ s 



x) 



P, 



(0) 



x)+\]'^^^,{x)\ + 0{al) (3.108) 



(3.109) 



where Se is defined in ( |2.3| ) and the coefficient functions C*-^^ and the a^^^ are given in 
appendix D. We have expressed the order corrections as sums of the NLO coefficient 



functions plus terms involving the a^^^, respectively a^'g. that will cancel with the 0{e) 
corrections to the lower order terms. The details of these cancellations will be dis- 
cussed in the next subsection in connection with the renormalization and factorization 
procedures. Only the class D of the ^I 'q has been calculated so that above equations 



.(1) 
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are determined only for the contributions with color factor NpT^Cp as indicated by 
the addition of an unknown contribution U that does not depend on Np. For this 
color class, as well as for the LO functions, we find perfect agreement with [j^,^. The 
complete set of functions C^^"* is given in appendix D. 

The color class F has been calculated here for the first time. It receives only 0{a1) 
contributions so that the result is regular. 



(3.110) 




Nc (Tq dx 



In the following subsection, the bare cross sections will be factorized in the MS- 
scheme and the coupling constant renormalization will be performed. 



3.3.3 MS Factorization and Renormalization 



As the cross section ( |3.110| ) does not contain singularities, it does not get factorized. 
The coupling constant renormalization is also trivial as this is the leading contribution 
to da^ and therefore the substitution (see, e.g. [63|) 



2^ 



27r 



2^ 2.^^UL 



(3.111) 



just introduces the scale ^r, 



1 



Ncctq dz 



' «s(/i|) 

27r 



(3.112) 



where now the variable z is used, just to be compatible with the notation in (|3.5| ). For 
the nonsinglet part of class D, the coupling constant renormalization actually serves to 



remove the 1/e singularity and the term proportional to a^l\. The LO functions C 



and ^a^^g appear in the subleading corrections of ( |3.107| ) with a common factor, which 



L,q 



is (apart from a {-) sign) the same as in ( p.lll| ). In the case of the pure singlet 
the renormalization only introduces a scale dependence in the strong coupling. The 
additional terms in (|3.108| ) are again proportional to the LO, but they are now removed 
by the factorization with the transition functions ( ^.4| ) and ( f^.8| ) according to ( |2.11| ). In 
the MS-scheme, the finite terms In {4tt) — 7e arising from the expansion of the factor 
are absorbed together with the singularity. Again just caring about terms proportional 
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to Nf in the 0{a 



2) 



1 dal^g _ asifi 



NcCTq dz 



1 

Ncao dz 



27C 



r) ^NS,{l)i 



4) 



27r 



2n 



1 r 
2 



1 



f^) In 



C 



ATS, (2) 
L,q 



c 



^^5,(2), 



(3.113) 
(3.114) 



where the logarithmic terms come from the expansion of the factors (/x^/ /i/j)"^, (/i^/ Q^Y, 
and {fi'^/Mjy, respectively, in powers of e according to 



l + e\n{z) + 0[e^) . 



(3.115) 



The complete set of longitudinal partonic cross sections is given in 0. There, the 
factorization and renormalization scales are identified. They will be disentangled here 
by the replacement 



a,(Mj) ^a,(/i^) 



\ «s(/ifl) /3o^^ 
27r 2 



(3.116) 



The perturbative expansion is performed in powers of ag/AiT in [0, whereas in this work 
the expansion parameter is as/27i. In our notation, the 0{al) coefficient functions thus 
get an extra factor of 1/4. The expressions are given explicitly in the appendix. The 
cross sections are defined by 



1 

NqCTo dz 



2tt 



{c 



NS,{2) 
L,q I 



+C'p ( 2 ln(l -z)~ \n{z) + - + z)\n 



1 

Nc<Jo dz 



1 dalg 
Ncdo dz 



+ 



NcCf ( - 



11 



+ NfTrCf 



In 



\f^R/ 



(3.117) 



27T 

+NfTrCf 



41n(;z) H 4:Z+ -z^ 

3 z 3 



In 



1 dal'^ ^ 1 dal^^ 



Ncctq dz Ncctq dz 
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(3.118) 
(3.119) 



1 da 



L,G 



Ncao dz 



2n 



2n 



{cg'„( 



+ 



Cp (^41n(z) -2 + 4- -2z 
+NcCf((-8 + 8-) ln(l 



NfTrCf + - 

^ 3z 3 



z - 



+ 8i)ln( 



Z 



38 461 



4 



+ 



H h 4^ z' 

3 3 ;z 3 

22 1 22 

J'z^l, 

8 1 8 

3^ ~ 3 



In 



Ml 



NcCf{-—- + — 
+NfTrCf 



In 



(3.120) 



in addition to eq. ( |3.112D . The convolutions with the C*^^-* (cf. ( |3.114| )) and the terms 
arising from the renormahzation (cf. (p.ll3| )) are evaluated and given as functions of 
z in above equations. When the scales are identified with Q^, as in the fits discussed 
below, the equations ( |3.112| ) and (|3.117| ) - ( |3.120| ) reduce to the generic form 



1 daY^^^ 



L,a 



NqCTq dz 



2ti 



-C 



label,(l) 
L,a 



(^) + 
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Ct:'^'\z) + 0({a,{Q')f). (3.121) 



The correct implementation of the coefficient functions in our computer program has 
been checked numerically by comparison of eq. (|3.8|) with the integral 



■ee,L 



= I dzz 



-'L,q 



' 2 



(3.122) 



3.4 Discussion 

After thus having reassured ourselves of the correctness of our numerical results, we 
will conclude with a short discussion of the NLO corrections to the longitudinal cross 
section. 

First, we are now in a position to assess the size of the K factor for the NLO. In 
the spirit of the definition of the NLO K factor as the ratio of the NLO result to the 
LO result, one might start to calculate the K factor for the differential cross section 
with the use of the NLO FF's in the numerator whereas the denominator is evaluated 
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Figure 12: The K factors for the longitudinal cross section, evaluated with our new charged 
hadron set (solid) and with the charged hadron set, as constructed in section 4.1 (dashed), 
a) shows the ratio of the cross sections evaluated consistently in NLO and LO, respectively. 
In plot b), the NLO FF's are used in both the numerator and the denominator and only the 
coefficient functions are changed. 



with the LO ones. However, as the two sets (LO and NLO, respectively) are fitted 
independently, they may differ significantly. In consequence, the K factor will depend 
strongly on the set of fragmentation functions used. This is demonstrated in Fig. |T2|a 
where the dashed curve shows the K factor obtained with our set of charged hadron 
FF's and where the solid curve has been calculated with our new set of charged hadron 
FF's. Because we used the longitudinal cross section measurement in our new fit, both 
the leading and the next-to-leading results describe the data well, which leads to a K 
factor close to one, whereas the K factor for the older set ranges between 1.45 at small, 
and 1.1 at large momentum fraction. It is thus advisable to keep the FF's at the same 
order in both the numerator and the denominator of K, to assess the size of the NLO 
corrections to the coefficient functions. This has been done in Fig. |K|b. The K factor 
is now insensitive to the set of FF's and is found to vary little, between 1.6 and 1.8. 
Note that the naive results of Fig. |T2| a would have been misleading with respect to 
the size of the NLO corrections. Fig. |l3b, on the other hand, confirms the conjecture 



prompted by Ree,L that the K factor would turn out to be large. 

Secondly, we show the individual contributions to the NLO longitudinal cross sec- 
tion, differential in x, in Fig. |1^. The results are obtained at the Z pole with our 
new set of fragmentation functions. The NS quark combination (solid) dominates for 
X > 0.2, whereas for smaller x the gluon (dashed) is more important. The PS quark 
combination (dash-dotted) and the contributions F (dotted) from the class CfTr only 
appear at NLO. For not too small x, the F competes with the negative PS (the ab- 
solute value is plotted). Whereas the relative weight of the pure singlet rises to a few 
percent at low the contribution of class CfTr is well below 1% everywhere. More- 
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Figure 13: The individual contributions to the NLO longitudinal cross section at the scale 
Q = mz, normalized to 1. The NS quark contribution (solid) dominates for x > 0.2, whereas 
for smaller x the gluon (dashed) is most important. Both the F (dotted) and the PS (das- 
dotted) quark contributions are absent at LO. At NLO, their impact is small and the F 
contribution competes with the negative PS (the absolute value is plotted) for not too small 

X. 



over, the relative importance of F decreases when moving away from the Z pole due 
to the decrease of the effective charge Q^{Q^), c.f. Fig. H. It can thus be neglected. 

The longitudinal polarized part contributes only a small fraction to the unpolarized 
(L + T) cross section, c.f. Fig. ^. With respect to measuring the gluon FF, however, 
it has the advantage that the gluon appears on an equal footing with the quarks in 
( p.38| ), whereas its contribution to ax is subleading. The gluonic contribution to ai is 
also numerically of the same order as that of the quarks, for not too large x, as can be 
seen in Fig. |13|. 

The longitudinal cross section, therefore, is a promising observable for a determi- 
nation of the gluon FF in e~^e~ annihilation. The error on the gluon FF has a huge 
impact when other processes are considered. This is due to large contributions from 
subprocesses with initial-state gluons that are usually associated with final-state glu- 
ons. Fig. |14| shows the ratio of the cross section, differential in pt that is due to 
gluon fragmentation in pp collisions at the TEVATRON and in photoproduction at 
HERA. In particular in pp scattering, the gluon FF dominates even for moderately 
large transverse momenta pt- For meaningful predictions of the IPP cross section in 
these processes the gluon must be well determined. This can be achieved by the use of 
the longitudinal cross section in the fit. 
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Figure 14: Fraction of the cross section that is due to gluon fragmentation (for the set |13| 
of FF's) for various processes and kinematics. The results for CDF and UA5 kinematics in 
pp cohisions at 1800, 630, and 200 GeV are shown by the sohd, dash-dotted, and dashed 
curves, respectively. The dotted curve gives the fraction for HI photoproduction kinematics. 
The drop at 3 GeV reflects the opening up of the additional charm channel. 
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4 Fragmentation Functions 



The investigation of fragmentation functions goes back to R.D. Field and R.P. Feynman 
who devised the first physical picture for the fragmentation process For a long 



time - from the early eighties to the early nineties - only LO fragmentation functions 
for charged pions, presented by R. Baier et al. pSl and updated by M. Anselmino 



et al. have been available. In retrospect, it is not easily understood why the 

progress in the field of parton density functions during that time-period has not been 
accompanied by similar advances in the field of fragmentation functions, despite the 
accumulation of a wealth of data. Whatever the reasons, the situation has changed 
for the better in the past few years by the presentation of numerous NLO sets of 



FF's. The first work was done by P. Chiappetta et al. ||6^ who extracted a variety 



of TT sets from e'^e annihilation, fixed-target, and collider data, supplemented by 
data for e^e~ annihilation obtained from a Monte Carlo simulation. We presented 



NLO sets for charged pions and charged kaons for the first time [^. Other, recent, 
sets are those for D*^ mesons obtained from e^e~ annihilation data by M. Cacciari 
and M. Greco and for photons by L. Bourhis et al. ||6^. Apart from these sets. 



and the ones by the author and collaborators G. Kramer and B.A. Kniehl pr3|JT^Jl5 
that are discussed below, all other sets we are aware of have been fitted to data from 
Monte Carlo simulations of e^e~ annihilation. Sets for the r] meson, for charged pions, 
and for charged and neutral kaons have been obtained by M. Greco and collaborators 
BP]. P. Nason and B.R. Webber have concentrated on the sum over charged hadrons. 



parametrizing the FF's at the Z pole, and also partly incorporated experimental data 
directly See also [^] for a review of recent work on fragmentation functions. 



Before discussing the new set of charged hadron FF's, as well as the sets pr3|JI^JT3 
that we published earlier, in subsections 4.1 through 4.4, some general aspects will be 
covered in this introductory section. Our general policy for both the improved charged 
hadron fit performed in this work and for the earlier fits is outlined in the following. 
Some exceptions that apply to the individual analyses are specified in the relevant 
subsections. 

As advocated earlier, we use only experimental data from e"'"e~ annihilation to fit 
the FF's, in order to preserve the most of the factorization theorem's predictive power. 
For its good statistics, data from LEPl are selected in all of the analyses presented. In 
some of the fits, data at 29 GeV are used in addition to the 91.2 GeV data to better 
describe the scaling violation. The benefits of this addition will be demonstrated in 
subsection 5.2. The hadrons are treated as massless so that this approach becomes 
unreliable for x Ivuhj \fQ^ and a lower limit around Xmin ~ 0.1 is enforced for the 
data used in the fits, as well as for subsequent applications of the FF's in section 5. 
Data above Xmax = 0.8 are also discarded for poor statistics and limited reliability 
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of the systematic errors in the experimental data. Another problem related to the 
data is that the phenomenology of fragmentation given here does not apply to hadrons 
originating from weak decays. Experimentally, such contributions can not be separated 
from the strongly produced hadrons at present. However, the corrections are of order 
rv?IQ'^ and amount only to about 0.1% at the Z pole They are consistently 

neglected in this work together with all other possible power-law corrections. 

The cross sections have been calculated to NLO only in the massless theory so 
all quarks are taken to be massless in this work. (An exception is the -D*^ analysis 
where the c and h quark masses are taken into account also by a modification of the 
factorization scheme.) We adopt the variable flavor number scheme for the treatment 
of the heavy-quark thresholds. The mass effects of the heavy quarks, c and 6, would be 
most pronounced in the threshold regions around 2mc, respectively 2mfe. Our massless 
approach gives reliable results for the c and h fragmentation only sufficiently above 
their respective threshold. 

Novel data from ALEPH 110,03 and OPAL W% at LEPl that discriminate between 



different flavors, enabled us to fit the FF's separately for each fiavor. In the analyses 
presented here, particles and antiparticles are not discriminated in the final state so 
that = for all g^. The superscript /i^ always refers to the sum of positively 
and negatively charged h in our notation. Some additional constraints from SU(3) 
symmetry etc. are used to reduce the number of free parameters, e.g. DJ[ = . 

The procedure is thus to start at a low scale /iq = a/2 GeV with Np = 3 quark 
flavors plus the gluon and to parametrize the FF of parton a to fragment into hadron 



h in the standard form 



D^(x,/i^) = Ar,x°"(l-xf'". (4.1) 

The three quark FF's and the gluon FF are then evolved to fic = 2 rric, where the 
charm fragmentation function, parametrized in the same form ( [4.1|) , is added. The set 
of 4 + 1 FF's is evolved up to fib = 2 nib, where the bottom quark is added and the now 
complete set is finally evolved to the scale Mff^ Although we are primarily interested 
in NLO sets of FF's, we simultaneously perform LO fits to check the perturbative 
stability. The evolution and the calculation of the cross sections are done consistently 
in NLO (LO) for the NLO (LO) sets. Details concerning the evolution are given in 
section 2.1 and the numerical techniques pertinent to the solution of the AP equations 
are discussed in appendix C. 

A low starting scale, as used in our analyses, has two advantages, (i) The relatively 
long evolution to scales where the FF's are compared to data ensures that those FF's 



^In [|T3 14 , the heavy quark masses were chosen as half that of their lowest bound states, i.e., 
2mc = nirj^ = 2.9788 GeV and 2mb = mx = 9.46037 GeV according to the 1994 tables of the Particle 
Data Group |7^. For the new fit (as well as in [^) the masses were set to round figures, nic — 1.5 GeV 
and rub = 4.5 GeV (respectively 5 GeV in p5[). 
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are physical, i.e. compatible with the evolution and do not retain gross features of 
the functional dependence chosen for the input distributions, (ii) With due caution, 
the sets of FF's so obtained may also be applied at low scales. This is desirable in 
particular in ep collisions and has been exploited in a number of phenomenological 
studies of HERA data, see section 5. An alternative approach is to use the measured 
FF's at the Z pole and evolve them backwards to lower energy scales. This avenue has 
been taken in . 



The evolved set of FF's is convoluted with the partonic cross sections according to 
( p.5|) to obtain the longitudinal or the unpolarized (i.e. the sum of the longitudinal and 
the transverse) cross sections. For the leptonic part of the cross section, the interference 
of the Z with the photon is fully taken into account and the electroweak parameters 
are taken to be [sin^ 6 w{mz)]Ms = 0.2319, Tz = 2.4959 GeV, and mz = 91.2 GeVQ 
for the numerical evaluation of the cross section. The strong coupling is consistently 
evaluated to two (one) loops for the NLO (LO) fits. The number of flavors varies from 
three at low scales to five at LEPl energies. This does not only affect the /3-function 
of QCD, but also the value of A, which exhibits discontinuities at the heavy quark 
thresholds. Details on the evaluation of A and in the variable flavor number scheme 
are given in appendix A. The scales are identified with the center-of-mass energy of 
the hadronic process, /i/j = Mj = \f(^. 

The QCD result is then integrated over the experimental bins and compared to the 
data. The goodness of the fit is measured by the per bin, where for a single data 
point, 

■ (4.2) 

The parameters of the fit, i.e. the Na-, ota and are optimized through repeated 
comparison of the evolved FF's to data, with the use of MINUIT . 



The functions so obtained merely give a phenomenological description of a small 
subset of e"'"e~ data - those that have been fitted. As a check they are compared to 
e^e~ data that have not been used in the fit. When these data are collected at a 
different CM energy, this serves also as a test of scaling violation. The purpose of the 
entire procedure is finally to make predictions with the help of these functions. To 
this end, it is desirable to have a parametrization of the FF's not only as functions of 
the momentum fraction x but also of the scale Mf. Such parametrizations have been 
obtained for the sets |ll3| , p!^ that are discussed in subsections 4.1 and 4.3, respectively^^. 



Basically, the functions are parametrized in the form ( |4.1| ), where the parameters are 



We used the actual parameters from the PDG ||7l|, namely [sin^ 6'vk(™z)]ms ~ 0.2315, Tz — 
2.491 GeV, and mz = 91.188 GeV, for our new analysis. 

A FORTRAN routine returning the values of the FF's of any parton to charged pions, charged 
kaons, charged hadrons or neutral kaons is available upon request from: binnewieOmail.desy.de 
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now polynomials in 



,4,3) 

V in(/i2/A2) ; 



Further details of the parametrizations are given in the relevant publications. 

The following sections will dwell on the particulars of the individual fits for charged 
pions, charged kaons, charged hadrons, neutral kaons, and D*^, respectively. In sub- 
section 4.1, the first sets of charged hadron FF's are discussed, i.e. those for charged 



pions and kaons [13] that can, with some assumptions on the residual contributions to 
the cross section, also be applied to indiscriminate charged hadron production. The 
treatment of charged hadron FF's will be improved by a novel fit to both unpolarized 
and longitudinal cross section data in the next subsection. Owing to the inclusion of 
the NLO longitudinal coefficient functions discussed in the previous section, this new 
set will feature a particularly well constrained gluon FF. Fragmentation functions for 
neutral kaons are presented in 4.3. The discussion of the D*^ mesons in sub- 



section 4.4 introduces some additional concepts in dealing with the heavy quarks, c 
and h. 



4.1 Charged Pions and Kaons 

The largest contributions to charged particle production in e^e~ annihilation originate 
from pions and kaons - quite naturally we chose those hadrons for our first NLO fits. 
In [0, the FF's were fitted to precise data from TPC at 29 GeV and successfully 



described e^e~ data taken at energies from 5.2 GeV up to 91.2 GeV. However, we had 
to assume that the valence-type quarks on one hand, and the sea-type quarks on the 
other hand, fragment in the same way. When novel data from ALEPH became available 
that discriminate between different quark fiavors, this limitation could be removed. We 
were able to obtain a set of FF's with well constrained individual flavors by fltting 
to both ALEPH and TPC data. Only this second set will be discussed, as it represents 
the updated and improved version of our earlier work. 

More speciflcally, we used data for charged pion and charged kaon production from 
ALEPH |Tl|| and TPC |p. In addition, the novel data from ALEPH were employed 
where they did not discriminate between the charged hadrons but instead distinguished 
three cases, namely fragmentation of (i) m, d, s quarks, (ii) h quarks and (iii) all flve 
quarks |]10|- That analysis was done on the basis of 40 pb~^ of luminosity, taken during 



1992 and 1993. Heavy flavor events were enhanced in the h quark sample with an 
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impact parameter tag and an event shape tag. The overall normalization error was 
estimated to amount to 1%. 

In order to determine the flavor differences of the charged pion and charged kaon 
FF's in our fit, the residual contributions to the charged hadron cross section must 
be modeled. Off resonances, they are mostly due to protons and antiprotons that are 
taken to behave like pions with a relative normalization of 

/(x) = 0.195 - 1.35(a;- 0.35)2 (^^_^^ 

that has been inspired by the differential cross sections for inclusive charged pion and 
proton/antiproton production |TT|. Surprisingly, this naive guess turns out to be quite 



good not only for the shape of the spectrum but also for the relative importance of 
the heavy flavor contributions. Very recently, the SLD collaboration presented data 
on flavor tagged IPP in Z decays [|73| that conflrm a close analogy between the heavy 
quark contributions to pions and protons for x above 0.1. Residual contributions from 
hyperons, charmed mesons etc. amount to less than 10%, presumably. Thus, the 
cross section is to a good approximation given by 



[^ + fi^)]-nr + -n-- (4-5) 



dx dx dx 



By the use of the novel ALEPH data together with eq. ( [1.5[ ) the b fragmentation 
could be determined. This represents a major improvement over our earlier analysis 
where we had to assume that the b fragments as the other sea quarks although it is 
known to exhibit a softer x-dependence due to its mass. Data from OPAL |jl2[ on 
Dq'(x), measured in three jet events, were used in the flt to address the notorious 
problem of the gluon FF. 

This leaves us with the u, d, s, and c quarks. The bulk of pions produced in strong 
interactions is due to u and d quarks, whereas the u and s quarks are most prominent 
in kaon production. Since we use data for tt^ and together with the new ALEPH 
data, this removes most of the residual freedom for shifting contributions between 
different flavors. Only the c quark can not be constrained well by this procedure as it 



is neither tagged in [10| nor dominant in either of the charged pion or charged kaon 
data. 

The functions were parametrized at the starting scale with the single, physically 
motivated constraint that valence-type quarks should fragment in the same way, i.e.. 



D:^=Df and (4.6) 
Dr=Dr. (4.7) 

As we flt data at different CM energies, A can also be determined via scaling violation 
and is kept as a free parameter. This gives 2x5x3 + 1 = 31 flt parameters. We 
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exclude data below x = 0.1 where the formalism is bound to fail, as well as data above 
X = 0.8 that are plagued by large systematic errors. 



set 


flavor 


N 


a 


/5 




u, d 
s 
c 
b 


1.15 (1.09) 
4.25 (3.48) 
3.99 (4.51) 
4.01 (3.60) 


-0.74 (-0.85) 
-0.77 (-1.03) 
-0.79 (-0.86) 
-1.03 (-1.13) 


1.43 (1.44) 
4.48 (3.90) 
4.78 (4.53) 
7.86 (7.12) 


G 


5.53 (6.57) 


-0.32 (-0.46) 


2.70 (3.01) 




u, s 

d 

c 

b 


0.31 (0.38) 
1.08 (1.12) 
0.81 (0.62) 
0.61 (0.73) 


-0.98 (-1.23) 
-0.82 (-0.92) 
-0.69 (-0.67) 
-0.88 (-0.80) 


0.97 (1.06) 
2.55 (2.85) 
2.98 (2.48) 
2.93 (2.83) 


G 


0.31 (0.37) 


-0.17 (-0.21) 


0.89 (3.07) 



Table 4: Parameters of the charged pion and charged kaon fragmentation functions obtained 
in the NLO (LO) fit. 



The values resulting from the fit are given in Table ^. Discussion of the results for A 
is postponed to section 5.1. Except for Dq^ , the parameters exhibit good perturbative 
stability. 



[GeV] 


7r± 


ref. 




ref. 




ref. 


91.2 


0.8 (0.9) 


|1 


1|* 


0.5 (0.8) 


[11|* 


0.4 (0.5) 


I 


rg* 


34.0 - 35.0 


1.3 (0.8) 


\74\ 


0.9 (1.0) 


1 


74 




0.4 (0.8) 




|75| 


29.0 


1.1 (1.3) 


1 


1* 


0.7 (0.8) 


i 




1.5 (1.4) 




78 1 


9.98 - 10.49 


1.5 (2.0) 


m 


1.1 (1.3) 


1 


76 


1 






5.2 


1.4 (0.5) 


|77| 


1.6 (2.1) 


1 


77 


1 







Table 5: P^r degree of freedom for the NLO (LO) fits of charged pions, charged kaons, 
and the sum over charged hadrons. The data used in the fits are marked by asterisk. 



The resulting from the fit, as well as those for a selection of other e^e date 



are given in Table BPl Good agreement with the data over a wide range of scales 



The values for the flavor separated sets of data are not listed, they are all less than 1 in NLO 
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is achieved. Most remarkably we find very satisfactory agreement even at low scales 
down to 5.2 GeV. The NLO (LO) cross sections are compared to e~^e~ data in Figs. |l^ 
and |TB|. Figs, ^a, b show the results for charged pion and charged kaon production, 
respectively. The sum over all charged hadrons is considered in Fig. The QCD 
predictions are in very good agreement with all the e~^e~ data, even for x considerably 
below Xmin = 0.1. Perturbative stability is evident also by the near coincidence of the 
NLO (solid) and LO (dashed) curves. As a further test, we checked the momentum 
sum rule; this has already been discussed in section 2.3. 











LO -— 




. NLO 





a) 



0.4 0.6 

X 



+1 

b 
-a 



b) 




Figure 15: Inclusive production of charged hadrons in e'^e annihilation. The plots show the 



production of a) pions and b) kaons. The data are from top downward DASP [77|, ARGUS 
||76l, TPC |], TASSO |7|, and ALEPH For better separation the data are multiphed 
by powers of 10. Solid (dashed) curves give the NLO (LO) results. 



This concludes the description of our fit [|l^. From the good agreement with data at 
lower CM energies we are confident that these FF's describe the e~^e~ annihilation data 
well and may sensibly be applied to other processes in section 5. The proton/ antiproton 
final state is not as prominent as in e~^e~ annihilation in other processes. We will thus 
approximate the subset of all charged hadrons just by charged pions and charged kaons 
in our applications to ep, pp and 77 scattering. 



with the exemption of the gluon data which could not be described well. 
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Figure 16: Inclusive production of charged hadrons in e'^e annihilation as in Fig. H but 



for all charged hadrons. The data are from top to bottom: MARK II [78|, CELLO [75|, 



AMY [79 1, and ALEPH |]lO| as well as u, d, s-enriched and 6-enriched samples fom ALEPH 



and as measured by OPAL y]. 

4.2 Charged Hadrons — New Fit 

In the following, the details of the new fit are presented. The approach is much the 
same as in the preceding subsection, i.e., we use the data from ALEPH in the fit0. 
The data from ALEPH enable a determination of FF's for individual flavors. In order 
to minimize the bias in the fit, we will not distinguish between pions and kaons but 
rather deal with the subset of all charged hadrons. 

In contrast to our earlier work, now also a charm-enriched sample is available in 
addition to the light-flavor and the bottom-enriched samples. We had to assume that 
the distributions are for pure samples in our earlier work. The values of the mixing 
matrix have been published in in the meantime. This enables us to construct pure 



samples by multiplying the data by the inverse of the mixing matrix. Fig. 17 shows 



Analyses of the flavor differences similar to that of ALEPH are currently beeing conducted 
by the DELPHI and OPAL collaborations. The longitudinal cross section that is also given in p5[ | 
has also been measured by OPAL |Q. The data from the ALEPH collaboration have been selected 
for their better statistics. 
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Figure 17: Flavor separated data from ALEPH |45] together with the results from the 
central fit. From top to bottom, the differential light quark, charm quark, and bottom quark 
distributions are shown, separated by factors of 10. 



the data resulting from this procedure which have been used for our fit. From top to 
bottom, the pure light (m, (i, s) quark, c, and b quark data are shown, separated by 
factors of 10. One point of the charm distribution (bin 0.75 to 0.8) turns negative in 
the demixing procedure. 

In addition to the unpolarized flavor tagged data we also use the data for the 
longitudinally polarized cross section that are shown in Fig. |TB. 



In connection 



with the NLO coefficient functions they help to constrain the gluon FF. 



data 






u, d, s 


0.14 


0.80 


c 


0.14 


0.75 


b 


0.20 


0.75 




0.08 


0.80 



Table 6: Cuts on the ALEPH data Q that are used in the fits. 



Only those bins were used in the fit, for which the correction factors in the experi- 
mental analysis were close to one. This restriction excluded some high-x bins. Lower 
limits were also enforced for theoretical reasons. For the longitudinal cross section, 
X > 0.08 was imposed, whereas the lower bound was stricter for the fiavor separated 
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Figure 18: The longitudinal cross section in e^e^ annihilation. The results from the fit are 
compared to the data from ALEPH Q. Both the NLO (sohd) and the LO (dotted) are 
in excellent agreement with the measurement, even at very low momentum fraction. The 



dashed curve shows the result from |13] at the scale niz- 



distributions. This can be justified because the small x behavior is dominated by the 
gluon in any case so that little information is lost by the exclusion of those points. The 
lower and upper bounds for the individual sets of data are listed in Table ^. 

The updated ALEPH data enable to fit light quark, c quark, and h quark FF's 
separately, without any further model assumptions. While we could reasonably assume 
that the u and the d quark behave similar in charged pion fragmentation, this does not 
hold for the sum over all charged hadrons. Therefore, we do not distinguish between 
the light quarks in this fit. We use the standard parametrization at the starting scale 
so that we have 12 fit parameters. 

With the extraction of as in ep collisions in mind, which will be discussed in sec- 
tion 5.2, we will not fit A but instead fix the parameter at certain values. As we do 
not want to fit A here, we do not include data at lower energies in our fit. We will 
present separate sets of FF's for A(^)'^-^° = 150, 200, 400 MeV. A^^)'^^'^ = 350 MeV 
corresponds to A^^^^'^^*^ = 229 MeV (see appendix A), which is very close to the value 
of 227 MeV we obtained in our earlier fit. This is our preferred choice, it will be referred 
to as the central value of A. As always, we will also present a LO fit for the central 
value. As the LO does not describe the scaling violation as well we will only present 
NLO sets for the other values of A. For the central choice of A, we also construct sets 
with /io = 1 GeV^ and /ig = 4 GeV^, in addition to our standard choice of /Xq = 2 GeV^. 
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order 


M, d, s 


c 


b 




average 


NLO 


0.91 


0.34 


1.63 


0.07 


0.78 


LO 


1.00 


0.35 


1.15 


0.09 


0.69 



Table 7: values for the central sets, fitted with A(4).^iO ^ 350 j^eV and = 2 GeV^. 



/x^ [GeV^] 


u, d, s 


c 


b 




average 


1 


1.25 


0.39 


1.63 


1.03 


1.06 


4 


1.06 


0.35 


1.99 


0.43 


0.98 



Table 8: values for the NLO fits with low and high starting scales, with A^"^)' 
350 MeV. 



A(4) [MeV] 


u, d, s 


c 


b 


o-L 


average 


150 


0.81 


0.33 


1.85 


0.18 


0.81 


200 


0.89 


0.32 


1.80 


0.24 


0.84 


250 


1.36 


0.40 


1.64 


0.64 


1.03 


300 


0.89 


0.34 


1.73 


0.16 


0.81 


400 


0.83 


0.34 


1.57 


0.11 


0.74 



Table 9: values for the NLO fits with the various values of A, and ^Ug = 2 GeV^. 



The obtained in the fit are given in Tables ^ through ^ for the central values of 
the parameters A^^^ and /iq, for the high and low starting scales, and for the various 
A values, respectively. Averaged over all data, we obtain values around Xd.o.f. — 0-8 
consistently for any choice of A or fiQ. 

The comparison of our central fit results with the ALEPH data have already 
been shown in Figs. ^ and |T8|. Both the NLO (solid) and the LO (dotted) describe 



the flavor tagged data well for x between 0.1 and 0.75. Our results are slightly low at 
higher x, whereas they overshoot the data for x below 0.1. Fig. |18| demonstrates very 
good agreement for the LO as well as NLO longitudinal cross sections, down to very 
low X. Because the NLO longitudinal coefficient functions were not yet available then, 
we combined the NLO FF's with the LO coefficient functions in [|l^. The result is also 
shown in Fig. as dashed curve. As the K factor for da/dx is large, c.f. Fig. |T^, the 
consistent use of the NLO coefficient functions would lead to considerably improved 
agreement even without a new fit. This was demonstrated in for our flrst set pT| . 
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order 


u, d, s 
N a p 


c 

N a P 


b 

N a P 


G 

N a P 


NLO 
LO 


2.05 -0.733 1.77 
1.79 -0.759 1.82 


4.22 -0.761 3.52 
4.52 -0.684 3.59 


2.39 -1.10 4.71 
1.81 -1.26 4.35 


5.47 -0.740 2.33 
4.30 -1.17 1.25 



Table 10: The parameters resulting from the central fits in NLO and LO. 



A{4) 

[MeV] 


u, d, s 
N a P 


c 

N a p 


b 

Nap 


G 

N a p 


150 
200 
250 
300 
400 


1.47 -0.950 1.88 
1.67 -0.970 1.96 

1.88 -0.693 1.81 
1.92 -0.842 1.85 

1.89 -0.766 1.64 


4.03 -0.766 3.73 
3.89 -0.817 3.64 
4.38 -0.740 3.78 
4.72 -0.738 3.75 
3.95 -0.802 3.43 


2.29 -1.07 4.72 
2.33 -1.09 4.72 
2.31 -1.07 4.68 
2.48 -1.07 4.77 
2.31 -1.12 4.66 


4.97 -0.935 2.33 
4.05 -0.876 2.20 
4.97 -0.992 2.28 
5.07 -0.729 2.30 
6.27 -0.561 2.31 



Table 11: The parameters resulting from the different choices of A in NLO. 



[GeV2] 


n, d, s 
N a P 


c 

N a P 


b 

Nap 


G 

N a P 


1 

4 


2.59 -0.638 1.70 
1.84 -0.851 1.98 


4.11 -0.839 3.47 
4.33 -0.708 3.54 


2.41 -1.11 4.67 
2.35 -1.05 4.64 


5.25 -0.710 2.32 
5.17 -0.683 2.24 



Table 12: The parameters resulting from the fits with low and high starting scales fiQ. 



Table |T0| lists the values of the parameters A^, a and P that result from the central 
fits. For further application, the parameters for the other values are also given in 
Tables |ll] and When comparing the parameters for the different choices it is 
apparent that the individual fragmentation functions are well constrained at NLO. 
In particular for the gluon, the results from the fits do not depend strongly on the 
choices of A or fiQ. This serves as prove that the gluon is indeed well constrained by 
the inclusion of the longitudinal data, at least to the point that it is not particularly 
strongly correlated with the other partons. As expected from their softer spectra, the 
values for the parameter P for the heavy c and b are larger than for the light u, d, s 
quarks. From around 1.8 for the light quarks, it increases to 3.7 for the charm and 4.7 
for the bottom quark. 
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a) b) 

Figure 19: The NLO fragmentation functions of the new fit. The spread for different 
values of the starting scale is shown at 10 GeV. a) The light quark (dashed) charm (dotted) 
and bottom (solid) FF's. The three curves for each parton are hardly distinguishable, b) 
The FF of the gluon. Going upward, the three solid curves correspond to values of //q of 1, 
2, and 4 GeV^, respectively. 



The sets for different choices of A represent different physics, namely different strong 
coupling strength so that the corresponding FF's are expected to differ. On the other 
hand, the starting scale of the evolution, fiQ, is an unphysical parameter and the results 
should not depend on it. The differences between fits for various values of /xq at the 
central value of A hence give a measure of how well the FF's are determined by the fit. 
Fig. p!9| a shows the results from the fits with /Iq = 1,2, and 4 GeV^ for the light quarks 
(dashed) the charm (dotted) and the bottom quark (solid). In Fig. |19|b, the gluon 
FF's of the same three sets are plotted as solid lines. While the quark FF's do not 
deviate much from each other, the gluon does exhibit a considerable spread at large x. 
Only at small x below 0.2, where it contributes dominantly to the longitudinal cross 
section (cf. Fig. |T3|), is it fixed by the fit. As the region of relatively small x is the 
most important, both in driving the evolution of the FF's and in the contribution to 
the cross section in hadronic or semi-hadronic scattering, this nevertheless represents 
a considerable improvement over our earlier analyses. 

Finally, we will compare the predictions of our new set of FF's to data at lower 
CM energies. Fig. ^ shows measurements of the unpolarized cross section taken by 
(from top to bottom) ALEPH |5[, MARK II |7|, and TASSO [||, at 91.2, 29 and 



14 GeV, respectively. The NLO results (solid) agree well with the data in the region 
of reliability between x = 0.1 and 0.75. The LO also gives satisfactory agreement. 
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Figure 20: The unpolarized cross section as evaluated with the new charged hadron FF's 
is compared to e~^e~ data at different CM energies. From top to bottom, the data are from 
ALEPH at 91.2 GeV ||, from MARK II at 29 GeV [^], and from TASSO at 14 GeV [|l|. 
Our NLO (LO) results are given by the solid (dotted) curves. For better separation, the data 
have been divided by powers of 10. 




X 



Figure 21: The ratio R of the new set of charged hadron FF's to those presented in [T^ at 
Mf = 10 GeV. Shown are the NLO results for the gluon (solid), the light quarks (dashed), 
the charm quark (dotted), and the bottom quark (dash-dotted). 
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We conclude that our new NLO set of charged hadron FF's features well constrained 
quark fragmentation functions. In addition to the bottom quark fragmentation, the 
charm quark fragmentation has now also been measured separately. Moreover, the 
gluon is better known than in our earlier sets, at least in the important small-x region. 
We will suggest how to use the improved knowledge of the charged hadron fragmen- 
tation functions to measure as and to determine the gluon content of the photon in 
section 5. A direct comparison of our new FF's with those of |jl3| for charged hadrons 
(pions, kaons, and protons) is shown in Fig. We plot the ratio of the NLO results 
at Mf = 10 GeV for our new charged hadron set over the one discussed in the previous 
subsection. The additional information on the purity of the samples influenced the 
shape of the bottom FF (dash-dotted) at large x. The new fit also profited from new 
information on the charm fragmentation not at our disposal in the earlier fit. Thus, 
both the marked increase of the charm FF (dotted) and the decrease of the bottom FF 
(dash-dotted) at large x are signals of improved data. The light quark fragmentation 
(dashed) is similar for both sets in the central region of x. The gluon, however, shows 
a significant increase at small x compared to our older set, induced by the improved 
treatment of the longitudinal cross section. 



4.3 Neutral Kaons 

Neutral particles do not allow for a precise determination of momentum in conven- 
tional experimental setups. The neutral kaons are studied by their decay products, 
experimentally. The main decay channel for the short-lived CP-eigenstate is into two 
charged pions and for a highly energetic the secondary vertex is located about 
3 cm from the primary one. The longlived eigenstate decays outside the driftchamber. 
Thus, it is the that is measured, and as it is an equal mixture of the strangeness 
eigenstates and 1^ this corresponds to measuring the average of those two states. 
In the following, we shall collectively use the symbol for the sum of and 
(or and 1^). For not too high momenta, when the localization of the secondary 
vertex becomes difficult, the signature is very clear, resulting in precise measurements 
even on a large background. Consequently, spectra of neutral kaons in pp collisions 
5B| , |5^ and in photoproduction P^JTP[ ] have been measured with small errors. 



In order to use these excellent data to test QCD, fragmentation functions for neutral 
kaons are required. This was our motivation for the work presented in [Q. We used 
e^e~ data from MARK II |^5|] and ALEPH ||8^. In the spirit of the earlier work on 



charged pions and kaons the standard parametrization (|4.1|) was used for each fiavor 
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individually. Only the valence-type quarks were identified, i.e. = Df . However, 
as no experimental information comparable to the one exploited in subsection 4.1 on 
flavor differences was available, the parameters were not well constrained by the fit. 
In particular, the gluon was strongly correlated with the other partons, which is a 
notorious problem in fits to e^e~ data. Inspired by isospin symmetry, we made the 
assumption that the gluons behave as in charged kaons, i.e. the respective FF's were 
identified at the starting scale, 

D^\x,ii,) = D^^{x,ii,). (4.8) 



The h quark is then fixed by the fake scaling effects that will be discussed in section 5.1. 
In essence, the change of shape from 29 to 91.2 GeV, after subtraction of the scaling 
violation effect which is dominated by the gluon, is attributed to the h quark. The 
other sea-type quarks are also rather strongly correlated with each other. A was kept 
fixed at the NLO (LO) value obtained in ||13| . 



10' 




0.0 0.2 0.4 0.6 0.8 1.0 

X 



Figure 22: Differential cross sections of inclusive + production at LO (dashed lines) 
and NLO (solid lines) as functions of x at ^/s = 91.2 and 29 GeV. The theoretical calculations 
are compared with the respective experimental data by ALEPH |86| and MARK II |85]. For 
better separation, the distributions at 29 GeV have been divided by 10. 
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The quality of the fit to the MARK II [|5l and ALEPH data is demonstrated 
in Fig. Both the NLO (solid) and the LO ( dashed) fit the data very well. The 
resulting from the fit are listed in Table O. Besides the ones used in the fit, data at 



lower CM energies also yield small values. 



[GeV] 


NLO 


LO 


ref. 


91.2 


0.5 


0.5 




86 


* 


35.0 


0.2 


0.2 




87 


1 


29.0 


0.4 


0.5 




85 


* 


10.49 


1.1 


1.4 




88 





Table 13: P^r degree of freedom for the NLO and LO fits of neutral kaons. The data 
used in the fits are marked by asterisk. 



set 


fiavor 


N 


a 


/3 




u 


0.53 (0.54) 


-0.57 (-0.77) 


1.87 (1.49) 




d, s 


1.45 (1.54) 


-0.62 (-0.72) 


3.84 (3.70) 




c 


1.70 (1.13) 


-0.51 (-0.70) 


3.76 (3.02) 




b 


0.47 (0.64) 


-0.66 (-0.63) 


1.49 (1.84) 




G 


0.31 (0.37) 


-0.17 (-0.21) 


0.89 (3.07) 



Table 14: The parameters obtained from the NLO (LO) fits for the neutral kaon. 



The parameters that result from the fit are given in Table |Tj. The perturbative 
stability is quite good for the differential cross sections, although there are some changes 
in the parameters of the sea-type quarks from LO to NLO. The overall agreement with 
e~^e~ annihilation data is good. The FF's will be applied to photoproduction and pp 
collisions in subsection 5.1. 



4.4 L>*± Mesons 



The HI [0 and ZEUS pO| collaborations at HERA presented new data on inclusive 
D*"^ spectra that extend up to px = 12 GeV. In order to obtain definite predictions for 
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£)*± production cross sections in the massless-charm scheme, one needs the PDF's of 
the charm quarks in the proton and the photon besides the FF's for D*"^ mesons. In |^ 



it was found that the charm component in the proton contributes only marginally to 
the cross section of inclusive D*^ photoproduction. This offers the possibility to use the 
new data from HERA together with NLO fragmentation functions for D*^ mesons to 
obtain information on the charm distribution in the photon. This was our motivation in 
[ITSj . This subsection contains the discussion of the sets of FF's presented there. Their 
application to the determination of the photon PDF's is postponed to subsection 5.3. 

Precise data on inclusive D*^ production in e^e~ annihilation have been presented 
by the ALEPH |Q and OPAL collaborations at LEPl, recently. The contributions 



from charm and bottom quarks are disentangled in those studies and light quarks are 
negligible in this process. With the use of these data, we obtained very well constrained 
sets of FF's. In set S we used the standard form ( ^.1| ) for both c and h quarks whereas 



we used the Peterson distribution [^] 



for the charm quark and the standard form for the bottom quark in a mixed set (M). 
The Peterson distribution is particularly suitable to describe FF's that peak at large 
X, so a good fit could not be obtained when applying it also to the relatively broad 
bottom FF. As the fragmentation of a light quark into a D*^ is very unlikely, the 
light quark FF's were set to zero at the starting scale. The gluon is generated only 
through the evolution as well. Therefore, the evolution starts at /io = 2mc = 3 GeV 
with four flavors. The h quark is added at 2mfe = 10 GeV. Values for A were taken 
from the charged pion and charged kaon analysis. We use the n-space evolution for set 
S whereas we found it more convenient to evolve set M in x-space. For a comparison 
of the two techniques see appendix C. 

So far, the formalism is identical to the one for describing the fragmentation into 
light hadrons based on the MS-factorization scheme. To account for the mass of the 
heavy quarks, we now adjust the factorization of the final-state coUinear singularities 
associated with the charm and bottom quarks in such a way that it matches the massive 
calculation. This is achieved by a change of the factorization scheme, i.e. we substitute 
P^alQ{z) ln(gVM|) ^ Palqiz) IniQ^M]) - dQa{z), in the coefficient functions where 
the dga are given by 

dQQiz) = -P^qIq{z) In + |-25(1 - ^) + 2 (- 



z 



+4 



ln(l - z) 



;i + 2)[l + 21n(l-^)] , (4.10) 
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iQg{z) = dQq{z) = dQQ{z) = , 



with Q = c,b and q = u,d,s. These functions have been extracted from 



(4.11) 
(4.12) 
m by 



B.A. Kniehl, G. Kramer and M. Spira in p9|, where a detailed discussion of the modified 
scheme is presented. The same substitutions must also be performed in any massless 
hard-scattering cross sections that will be convoluted with these FF's. When we study 
D*^ production in photoproduction in the next section, this has of course been taken 
into account appropriately. 



set 


flavor 


N 


a 






LO S 


c 


483 


8.04 


3.08 






b 


204 


2.96 


6.05 




NLO S 


c 


496 


8.84 


2.90 






b 


173 


2.81 


6.20 




LO M 


c 


0.214 






0.0926 




b 


197 


3.04 


5.92 




NLO M 


c 


0.161 






0.0674 




b 


194 


2.97 


6.22 





Table 15: The fit parameters for the charm and bottom quark FF's of sets S and M at LO 
and NLO. The other FF's are set to zero at the starting scale. 



The values of A^, a, /3, and ec that result from the LO and NLO fits to these 
data are displayed in Table |TB|. The parameter ec appears only in the charm quark 
FF's of set M. Bottom fragmentation is small at low energies and one may then use 
ec as the single parameter in phenomenological analyses of D*^. This parameter thus 
supersedes the value by Chrin |^ that was obtained from older data with different 
methods, ec depends on the scheme of factorization and on the way perturbative and 
non-perturbative effects are handled in the fragmentation functions. The parameters 
of the bottom quark FF's of sets S and M are very similar, i.e., the data fix the charm 
and bottom quark FF's independently of each other. 

The results of our fits S and M are compared to the ALEPH data |Q in Figs. p3| a 
and b, respectively. Figs. |2^ and b shows the analogous comparison for the OPAL 
data Except for very small x, the NLO (solid) and LO (dotted) curves are very 
similar. This is also true for the distributions at the starting scale, as may be seen by 
comparing the corresponding LO and NLO parameters in Table |T^. Only ec of set M 
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changes noticeably from LO to NLO. For a; — 0, we obtain large differences between 
LO and NLO, indicating that the perturbative treatment ceases to be valid. In this 
limit, also the massless approximation is not good anymore as pointed out already. 
Our results are meaningful only for x > Xcut = 0.1 and the first bin of the OPAL data 
is excluded from the analysis. 




Figure 23: Inclusive D*^ production in e^e~ annihilation. Our predictions for a) set S 
and b) set M are compared with data from ALEPH |9C]. The three sets of curves and data 
correspond to the Z ^ cc and Z ^ bb samples as well as their sum. The NLO (LO) results 
are plotted as solid (dotted) lines. 




Figure 24: As Fig. ||, but for the OPAL Q data. 



The Figs. P^a, b and P4|a, b show the normalized cross section ( p.5|) , which also 
includes the gluon contributions in NLO. These are attributed to the charm and bottom 
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quark curves, respectively, in proportion to their charge weight. The contributions from 
a = u,u, d, d, s, s are also included for the curves that give the total D*^ cross section. 
Numerically, their contribution is less than 2% of the integrated cross section. 

The values obtained with sets S and M at LO and NLO for the ALEPH and 



OPAL data, as well as their combination, are listed in Table |16[ Along with the total 
values, we also give those for the individual subsets of data. Set S leads to better 
values than set M, which is not surprising if one recalls that set S has six degrees of 
freedom while set M has only five. The ^ire smaller for the OPAL data than for the 
ALEPH data, in general. 



set 


total 


ALEPH 




OPAL 








sum 


c 


b 


sum 


c 


b 


LO S 


0.81 


0.73 


1.29 


0.65 


0.48 


1.02 


0.82 


NLO S 


0.79 


0.78 


1.31 


0.63 


0.42 


1.08 


0.63 


LO M 


1.03 


1.05 


2.11 


0.77 


0.47 


1.29 


0.75 


NLO M 


0.95 


1.10 


2.05 


0.74 


0.40 


1.11 


0.47 



Table 16: The P^r degree of freedom obtained in the LO and NLO fits S and M to the 
ALEPH @ and OPAL [|| data. The first bin of the OPAL data has been excluded from 
the fit. 



The integrals of the charm and bottom quark EE's into D*^ mesons over x give the 
branching fractions of these transitions. Eor the reasons given above, we restrict our 
considerations to the region Xcut < x < 1, with Xcut = 0.1, and define 



dxDnix, M 



Xcut 



f' 



(4.13) 



where Q = c,b. Experimentally, the contribution from the omitted region < x < 
Xcut is close to zero with a large error. When evaluating ([4.13|) with our sets S and 



M, we find that the values depend only weakly on the scale Mf. The theoretical 
results for the branching fractions are in good agreement with the values measured 
experimentally at LEPl. Eor example in NLO, the sets S and M yield 1.12 and 0.93 
respectively for the ratio of the charm over the bottom branching fractions. This agrees 
well with the experimental result of Bc{mz) / Bh{mz) = 1.03 ± 0.15 [|91|. In a more 
standard approach, the EE's would have been normalized to the branching fractions 
and eq. ( |4.13| ) would have reproduced these experimental results by construction. In our 



approach the agreement is nontrivial and indicates the soundness of the fit procedure. 
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The first moment of tlie FF's, normalized to tlie brandling fraction, gives tlie mean 
momentum fraction, 



(^)q(M/) 



dx xDq{x, M 



f' 



(4.14) 



where Q = c,b. The values of {x)Q{Mf) for Q = c,b evaluated with sets S and M in 
LO and NLO at Mf = 2m,Q, mz are collected in Table 0. The differences between sets 
S and M and between LO and NLO are small. The effect of the evolution, however, 
is significant, e.g., {x)c{Mf) decreases from about 0.7 at the effective starting scale 
to about 0.5 at mz- When these results are compared to the experimental numbers 
reported by ALEPH Q and OPAL m. 



, , , , +0.010 
(x)c(mz) = 0.495_^^^^± 0.007 

, , , , +0.008 
(x),(mz) = 0.515_^^^^^± 0.010 



(ALEPH) 
(OPAL) , 



we find good agreement for {x)c{fnz)- 



(4.15) 
(4.16) 



set 


{x)c{2mc) 


{x)c{mz) 


{x)b{2mb) 


{x)b{mz) 


LO S 


0.689 


0.518 


0.364 


0.320 


NLO S 


0.716 


0.528 


0.351 


0.309 


LO M 


0.637 


0.486 


0.372 


0.327 


NLO M 


0.664 


0.495 


0.359 


0.315 



Table 17: Mean momentum fractions of D*^ mesons produced in charm and bottom quark 
fragmentation. Eq. (4.14) is evaluated at the respective starting scales and at Mf = mz with 
sets S and M in both LO and NLO. 



As a final check we compare our FF's to e~^e data at lower CM energies. For the 

[g at = 10.6 GeV, from HRS |9| at 



comparison, we selected data from ARGUS 



^/s = 29 GeV, and from TASSO |9| at = 34.2 GeV. The TASSO collaboration 
measured two decay channels of the meson, namely, K^tt^tt'^tt" (TASSO 1) 

and D° ^ K-TT+ (TASSO 2). Fig. || shows the ARGUS, HRS and TASSO data on 
the cross section Q'^da/dx of + e~ D*'^ + X, together with the respective LO and 
NLO predictions based on set S0. The theoretical results are calculated with Np = 5 

-'^^For consistency, we corrected the HRS and TASSO data by updating the values of relevant branch- 
ing ratios to the new values from the 1994 tables of the Particle Data Group [Q, which were also 
used by OPAL The ALEPH data used in the fit have also been corrected accordingly. 
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quark flavors (except for tlie case of ARGUS, wliere Np = 4 is used) and also the 
contributions due to gluon and light-flavor fragmentation are consistently included as 
in Figs. ^ and For completeness, the OPAL data |^T| on the dimensionless cross 
section {l/(7ha_^)d(7/dx are also included in the plot. The agreement between the QCD 
predictions and the data is quite satisfactory. This ensures that the data indeed exhibit 
the scaling violation predicted on the basis of the LO and NLO evolution equations for 
our FF's. In fact, the change in the shape of the differential cross section with is 
mainly due to the bottom quark. The ARGUS data are taken off the resonance, at 
■\/s = 10.6 GeV, where the bottom quark is not yet active. The inclusion of the bottom 
quark contribution leads to a softening of the distribution, as may be seen in the case 
of the HRS data H at = 29 GeV and the TASSO data [||] at = 34.2 GeV. 
The evolution from 29 GeV to 34.2 GeV has no discernible effect. Going to the Z 
pole at y/s = 91.2 GeV further increases the relative importance of the bottom quark 



fragmentation and leads to an even softer spectrum. Except for the problematic region 
of small X, the NLO and LO predictions are very close to each other. 




Figure 25: Data for D*^ production in e^e~ annihilation at a number of CM energies are 
compared to the results obtained with the standard set of fragmentation functions. Solid 
(dotted) curves show the NLO S (LO S) results. For separation, the data have been rescaled 
by powers of 10. The topmost data (circles) are from OPAL [^]. The point at x = 0.05 is 
omitted and these data are displayed in the form l/dhaddf/dx. Going down on the plot and 



in energy, the other data are TASSO 1 (open triangles) and TASSO 2 (full triangles) |Q at 
34.2 GeV, HRS p| (squares) at 29 GeV and ARGUS [M (asterisk) at 10.6 GeV. 



The successful comparisons in Fig. E5| make us confident that our D*^ FF's, al- 



though constructed at y/Q^ = mz, also lead to useful descriptions of D*"^ fragmentation 
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at other scales. In the subsection 5.3, we shall exploit this property together with the 
universality of fragmentation to make predictions for inclusive D*"^ photoproduction 
at HERA. 

Recently, Martin et al. have presented an alternative method to account for massive 
quarks in parton density functions by modifying the evolution kernels as well as the co- 
efficient functions Their approach is well defined also in the problematic threshold 
regions and is therefore theoretically very appealing. An analogous treatment could be 
applied to the timelike regime and could represent a considerable improvement upon 
the analysis presented here. 
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5 Applications 



This section is dedicated to the apphcation of the fragmentation functions that have 
been extracted from e~^e~ data to a variety of processes in order to test QCD. First, 
universahty will be tested in a range of processes from ep scattering over pp scattering 
to 77 collisions. Then, we will investigate scaling violation, not only demonstrating 
its presence but actually determining Aqcd in the ir^/K^ fit. We will also discuss 
a consistent method to extract from p^ spectra of inclusive particle production in 
DIS and pp collisions. Finally, in section 5.3, we will look into the structure of the 
photon. Both its gluon and its charm content can be constrained with rapidity spectra 
for inclusive hadron photopro duct ion at HERA. 



5.1 Testing Universality 

As a consequence of the factorization theorem, the FF's extracted from e^e~ data in the 
previous section can be applied to predict IPP in any process. In this section, we will 
demonstrate this for ep collisions, pp collisions and 77 scattering. The photoproduction 
process will be discussed in some detail as it will also be used to study the PDF's of 
the photon in subsection 5.3. The discussion of pp and 77 collisions will be limited to 
the comparison of data with our predictions. 

In electron-proton collisions one usually considers two complementary scenarios: 
The limites of small virtuality of the photon - photoproduction - and of large virtuality 
- deep inelastic scattering. The former contributes the bulk of the ep cross section so 
that the experimental measurements benefit from high statistics. DIS, on the other 
hand, has the merit to probe the proton on the smallest scale in position space, giving 
the most detailed information about its partonic content. The intermediate region 
has not received much attention in the past but is presently being investigated in a 
number of new analyses |9^j9^. For DIS the situation is somewhat similar insofar as 
Pt spectra of inclusive particle production have not been analyzed to NLO yet, this is 
being changed at the moment ||10CI|| . A NLO comparison of the xp spectra with data 
from HI [ |101| has been carried out successfully by D. Graudenz ||102| | with the use of 



our set of charged particle FF's [|T3[- 

The Pt spectra for charged particles have been measured by the HI ||103|| and the 



ZEUS [104] collaborations. Very recently, the IPP process has been measured in DIS 



for the K'^ and D* final states, also [ p.05| , |106| . In this work, however, the DIS process is 
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Figure 26: Schematic diagram for IPP in ep scattering. 



disregarded in favor of the photoproduction process, where excellent data on inclusive 
particle production have been available for some years. 

Before comparing the NLO results on production with the data, the relevant 
kinematics are defined. The process can be depicted by the schematic diagram in 
Fig. |2^, where the virtuality is = —q^ ~ and the kinematics of the final state 
hadron are determined by its transverse momentum pt and its rapidity y. The CM 
rapidity, ycm, is related to yiab by 



1, (E^ 



(5.1) 



In the massless approach, the rapidity is equal to the pseudorapidity t] = — ln(tan(6'/2)). 
As usual, the photon spectrum of the electron/positron is approximated by the Weiz- 
sacker- Williams formula [p.07]| : 



a 



Mz) = - 




+ 2mp 



2 

max 



Q 



2 

min , 



(5.2) 



where z = E^/E^. Q^n = ^2^7(1 - z) and Q'i^^ = 0.01 GeV^ (0.02 GeV^) for HI 
(ZEUS) tagged events and Q^ax = 4 GeV^ for untagged events (both HI and ZEUS). 

It is well known that photoproduction proceeds via two distinct mechanisms. The 
photon can interact either directly with the partons originating from the proton (direct 
photoproduction) or via its quark and gluon content (resolved photoproduction). The 
cross section for the resolved process is given by 



E, 



d^a{-fp ^ h + X) 



E 

a.b.c 



Xu 



Ml)D^^{xu,M}) 
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X ■ 



+ ^4!^Kat^c{s, V, w- fil, Ml Ml Mj)j , (5.3) 
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where da'^f^^^/dv are the LO hard-scattering cross sections, and where v = 1+i/s and 
w = —u/{s + t) are functions of the Mandelstam variables at the parton leveL Those 
are defined as s = (p^ + Pb)^, ^ = {pa — PcY, and u = {pb — PcY- The parton momenta 
are related to the photon, proton, and hadron momenta by pa = x^q, pb = XpPp, and 
Pc = Ph/xh- The indices a, 6, and c run over the gluon and Np flavors of quarks and 
antiquarks. The Kab^c functions may be found in Ref. [|108|| for M^ = Mp. The G 



7 



and Gl are the PDF's of the photon and the proton, respectively. 

The NLO cross section of direct photoproduction emerges from eq. ( ^.3|) by sub- 
stituting G'^i^x^, M^) = 5(1 — x-y), replacing daab^c/dv and Kab^c with da^b-^c/dv and 
K^b^c, respectively, and omitting the sum over a. Analytic expressions for the K^b^c 
functions are listed in Ref. | [L09|| . 

When comparing to data, ( |5.3| ) is integrated either over rapidity or over the trans- 
verse momentum p^, leading to singly differential pt or rapidity spectra, respectively. 
The integration boundaries are fixed by the experimental conditions, details are given 
in our respective publications. Except where stated otherwise, we employ the photon 
PDF's of GRV flTol. In our recent work on D*^ production we used the set CTEQ4M 



111|| for the proton PDF's, whereas we used CTEQ3M ||112|| in our earlier analyses. 

Very precise data on IPP in photoproduction have been measured by the two HERA 
collaborations, ZEUS and HI. The first data to be published have of course been 
those for indiscriminate charged hadron production, as they have the best statistics 



1 8, 2 1,113 1 . Since then, more exclusive analyses have followed suit; for neutral kaons 



m 



rg and for D*^ in [0,^. 

The best data on charged hadron production to date are preliminary ones from HI 



113|| . The theoretical predictions possess an inherent ambiguity due to the unphysical 
factorization and renormalization scales. Before comparing to the high precision data, 
we will try to reduce this ambiguity by additional conceptual input. The cross section, 
when calculated to all orders, must not depend on the unphysical scales. The ambiguity 
introduced by the scales in this sense reflects our ignorance about unknown higher 
order terms in the perturbative expansion. With the prejudice that higher orders in 
the expansion should give only small corrections relative to the NLO, and disregarding 
for the moment that the expansion in ag represents only an asymptotic expansion, one 
arrives at the conjecture that the unknown corrections vanish for those scales where 
already the first nonleading corrections (NLO) disappear. This is known as principle 
of fastest apparent convergence. As a first approximation all the scales are identified 
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Figure 27: Scale dependence of the photoproduction cross section |114]. In departure from 
the usual choice, the proton PDF's of GRV |115| ] are employed in this plot. We use the 
charged hadron FF's of |13]. 
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Figure 28: Photoproduction for ^ = 0.6 compared to preliminary HI data The curve 



shows the NLO result for ^ = 0.6 ]116|]. The PDF's and FF's are the same as in Fig. 27 
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and parametrized by 



Mf = M^ = Mp 



(5.4) 

< 2.5 is plotted as 



The cross section integrated over > 4 GeV and —1.5 < yiab 
a function of the single parameter ^ in Fig. |2^. The scale dependence is markedly 



reduced from LO (dashed) to NLO (solid). Following the principle of fastest apparent 
convergence, we choose ^ = 0.6 for the scales in the comparison with the pt spectrum. 
Our curve is plotted together with the HI preliminary data in Fig. The agreement 



with the data turns out to be perfect. This will be exploited further to constrain the 
gluon density of the photon in section 5.3. However, we want to point out that the 
principle of fastest apparent convergence is just one of several methods for choosing 
a scale. It does not necessarily give the best choice of scales with respect to the 
convergence behavior of unknown higher orders. 
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Figure 29: Comparison of the preliminary HI tagged data on neutral Kaon photoproduction 
with the NLO prediction [|114| . Solid, dashed, and dash-dotted lines correspond to the NLO 
results with = 1/2, 1, and 2, respectively. 



Our NLO predictions for neutral Kaon production in 7p collisions are compared 
to data from HI (no longer preliminary) in Fig. The solid line gives the 
result for ^ = 1 whereas the dashed (dot-dashed) curves correspond to the choices 
^ = 1/2 (^ = 2), respectively. It turns out that the central value agrees best with 
the data. Within the theoretical uncertainty simulated by the three choices of scales. 
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good agreement is found down to very low scales oi px = 1.5 GeV where we would not 
expect our results to be reliable. 




Figure 30: The LO and NLO predictions of inclusive D*^ photoproduction in ep collisions 
based on sets S and M are compared with the ZEUS data We consider a) the px 

distribution da/dpT integrated over — 1.5 < yiab < 1 and 115 GeV <W< 280 GeV and b) 
the 2/iab distribution da/dy\ah integrated over 3 GeV < < 12 GeV and 115 GeV <W< 
280 GeV. 

The results of our recent work on D*^ mesons are shown in Fig. ^0|a, b. They are 



in good agreement with the px spectrum measured by ZEUS down to p-r = 3 GeV. 
Below this scale, the theoretical prediction is zero as the FF's start their evolution at 
2mc, effectively. Integration over pt yields the r] spectrum of Fig. |30[b. Within the 



large errors, we again agree with the data. We want to point out that the results of 
both sets, S and M, are displayed in Fig. and b. However, the LO curves of sets 
S (dotted) and M (dashed) almost coincide with each other, as do the NLO curves 
(dash-dotted and solid, respectively). In fact, they are absolutely indistinguishable in 
pO|a. This serves to prove that our results do not depend on the functional form of the 
start parametrizations. 

We will now turn to pp scattering. It is of particular interest to apply our FF's 
to hadronic interactions. As demonstrated in section 3.4, most of the pp cross section 
comes from the fragmentation of the gluon. This process may thus be used to test the 
gluon FF which is only weakly constrained by our older fits to e^e~ data. 



The IPP cross section has been calculated to NLO in |117| for the first time. With 



the use of our older set |31|, good agreement was found with the data from UAl- 
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Figure 31: Inclusive particle production in pp scattering. Our NLO QCD results [114, llC] 
are compared to data from UAl-MIMI at 630 GeV. The dashed, solid, and dash-dotted 
curves correspond to scale choices of ^ = 1/2, 1, and 2, respectively, a) Charged hadron 



production, calculated with our set |13]. b) Neutral kaon production with the FF's of |14]. 



MIMI,CDF, and CDHW [S4,rT|]. The NLO predictions that are tested against data 



here, have been calculated in [|116|| for the charged particle spectra and in [jT4| for 
the neutral kaon spectra. Fig. pl] a shows inclusive charged particle production in pp 
collisions. The NLO results with our set |13| are compared to data from UAl-MIMI 
8^ . The dashed, sohd, and dash-dotted curves correspond to scales ^ = 1/2, 1, and 
2, respectively, where ^ is defined in (p.4|). The lower choice of scales gives satisfactory 



agreement with the data over the whole pr range. The dotted curve will be explained 
in the next subsection. 



A similar comparison is presented for neutral kaon production in Fig. pTlb. The data 
are from UAl-MIMI also. Due to the more exclusive final state, the spectrum 
does not extend as far as for the charged hadrons. When comparing with the NLO 
result for our set of neutral kaon FF's we find again good agreement, this time for 
the central choice of scales. Our results are slightly high only for pjn below 1.5 GeV, 
where we do not expect our formalism to be valid. 

The last process that will be, briefiy, discussed is charged hadron production in 
77 collisions. So far, there exist only few data. However, with the prospect of future 
measurements at LEP2 and at the NLC, this situation is bound to improve. Here we 
present the results of ||23| that draw on the calculations presented in [|109|| . The process 



has been measured by the TASSO collaboration at PETRA | |119[ at a CM energy of 
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Figure 32: Inclusive diarged-liadron production in single-tagged 77 collisions. The 
MARK II data |l2|] = 29 GeV) are compared to the NLO calculations of da/dpT 

for scale choices ^ = 1/2, 1, and 2 p3[| . 



33.1 GeV, and by MARK II at PEP [ng at a CM energy of 29 GeV. We compare the 
MARK II data with our results in Fig. |32|. As usual, three curves for different choices 
of scales are plotted. The QCD result is too low, especially at high px- The same 
is observed for the less precise TASSO data (not shown). This feature has also been 
observed by other groups, e.g. |[109|| , but is not understood, yet. 

We now move on to make predictions for the future colliders LEP2 and NLC. Of 
the various experimental scenarios for the latter we will only consider the so-called 
LASER mode where laser light is Compton-backscattered off the e"*" and e~ beams. In 
Fig. ^ we display the relative importance of the various contributions to the spectra of 
single tagged charged particle production at LEP2 and at the NLC, respectively. The 
curves give the NLO predictions with our set [O at central scale ^ = 1. For LEP2, 
the direct -direct contribution is dominant for > 3 GeV [|2^ . As a consequence, the 
impact of the gluon, both in the PDF's and in the FF's, on the cross section is small 
as seen in Fig. p3|a. Switching off either gives the dashed, respectively dot-dashed lines 
that rapidly approach 1 with increasing px- The fraction due to the kaon is also plotted 
(solid), it increases slightly with p^, to about 0.4 at pr = 30 GeV. 

For the LASER mode of the NLC, on the other hand, the resolved-resolved contri- 
bution dominates up to pt = 20 GeV. This process proceeds mainly via gluon-gluon 
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Figure 33: Inclusive charged-hadron production in single-tagged 77 collisions, a) The 
influence of the gluon PDF of the photon, the gluon FF, and the charged-kaon final states. 
Shown are the fractions that remain if the gluon is switched off in the photon PDF's (dashed 
line) or in the FF's (dot-dashed line) as well as the fraction due to charged-kaon production 
(solid line), b) Influence of the gluon PDF of the photon and the gluon FF on the total result 
for the laser spectrum. Shown are the fractions that remain if the gluon is switched off in 
the photon PDF's (dashed line) or in the FF's (dot-dashed line) as well as the ratio of the 
calculation with the GS photon PDF's to that with the GRV set (solid line). 



fusion. Consequently, the gluon in the photon is much more important for the NLC 
The dashed line in Fig. p3| b shows the fraction of the total cross section that remains 
when the gluon in the photon is switched off. The solid line gives the ratio of the NLO 
results obtained with the GS photon PDF's | 121 | to those obtained with the GRV set. 
Both curves reveal a significant sensitivity to the photonic PDF's, in particular of the 
gluon. One may hope to exploit it to improve our knowledge of the gluon density in the 
photon. Unfortunately, this sensitivity is accompanied by a sensitivity of comparable 
size to the gluon FF, which constitutes a large part of the cross section (dash-dotted). 
Therefore, the latter has to be known well, for a meaningful measurement of the gluon 
in the photon. 
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5.2 Scaling Violation 



The high-precision data acquired at LEPl are, in union with lower energy data, an 
excellent starting point for the determination of the strong coupling. In fact, the as 
measurements obtained at LEPl are among the best to date. Of the numerous an- 
alytical methods employed, like event-shape variables or three jet over two jet rates, 
we will deal only with scaling violation of fragmentation functions. A dedicated anal- 
ysis along this line has been conducted by the ALEPH collaboration [^[^ that yields 
Us = 0.126 ± 0.009. We will take a more qualitative view and restrict ourselves to the 
extraction of Aqcd without a quantitative assessment of the theoretical and experi- 
mental errors. 

One common problem in the determination of the strong coupling from IPP in e~^e~ 
annihilation is fake scaling violation from the b quark contributions. The data with 
the best statistics are from LEPl so that one usually uses those together with lower 
energy data in the analyses. At the Z pole, the b is strongly enhanced, whereas it 
is negligible at lower energies. Due to its mass, the spectrum of b quark fragmentation 
is much softer than that of the light quarks so that an enhancement softens the IPP 
spectrum, mimicking scaling violation. Fortunately, this contribution to the change of 
shape is under good control in our set [|T3[ because information on fragmentation of 



the b quark was used in the fit, as detailed in section 4.1. 

We employed the scaling violation effect in our fit to data from 29 and 91.2 GeV, 
to determined Aqqj-, as one of 31 parameters. The values obtained from the fit were 
0.227 and 0.108 in NLO and LO, respectively. This translates into values of the strong 
coupling at the Z pole of Q;s(m|) = 0.118(0.122) at NLO (LO). The agreement of both 
the LO and NLO results with the value of 0.120 ± 0.008, extracted from a global fit to 
the observables measured at LEPl ||123|] , is striking. 



In our new fit we chose an alternative approach. Instead of aiming to determine 
from e+e~ annihilation data, we are interested in the DIS and photoproduction data. 
Compared to e+e~ , the ep process has the advantage that data at different scales 
are collected by the same detector in one measurement. With other words, the pr 
spectra are sensitive to Aqcd and can be used to measure a^. This is demonstrated 
in Fig. pria for pp scattering. The dotted curve results from the NLO calculation with 



^ = 1 when the evolution is switched off in the fragmentation functions. This leads 
to a significantly harder spectrum and is clearly inconsistent with the measurement. 
Thus the shape of the curve will be affected by a variation of A. In photoproduction, 
a change of 50 MeV in A*^^^ results in a 10% change of the NLO cross section at high 
Pt- 

For e^e~ annihilation, on the other hand, data at different experiments have to be 
combined to measure the scaling violation. This introduces large relative normalization 



'A similar analysis has also been conducted by the DELPHI collaboration |122|. 
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errors. They impair the potential benefit especially of data at low which would give 
a large scaling violation effect when combined with LEPl data. 

For a consistent fit of A to DIS or photoproduction data, A has to be varied not 
only in the coefficient functions, but also in the PDF's and in the FF's. PDF sets 
with variable A have been presented for the proton by A. Vogt ||124||, some time ago. 



No corresponding sets exist for the photon PDF's or for the FF's. The former are so 
poorly known that it does not make sense at the moment to construct a variable A set. 
The latter are presented in our new fit, for the first time. In the same manner as in 
121, we kept A^^) fixed at values of 150, 200, 250, 300, 350 and 400 MeV, respectively. 



in our NLO analysis. 
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Figure 34: Scaling Violation in e annihilation. The plot shows the ratio R of the 
differential cross sections for inclusive charged hadron production at 91.2 GeV to that at 
29 GeV. Data from ALEPH g|] and MARK II |7|l are compared to the predictions of 
the new /i^ fragmentation functions. The solid (dotted) lines correspond to the NLO (LO) 
results. 

After completion of the fit, we checked that the results describe the scaling violation 
in e^e~ annihilation correctly. The comparison of our central set (A^^) = 350 MeV) to 
data from ALEPH ^ and MARK II J7| at 91.2 and 29 GeV, respectively, is shown in 
Fig. 0. Both the NLO (solid) and LO (dotted) sets of FF's lead to good agreement with 
the scaling violation measurements. The same holds true for the sets with various A, as 
demonstrated in Fig. The results for the individual sets exhibit a steady decrease for 
successively higher values of A, also for the sets with A*^^) = 200, 250, 300 MeV, which 
are not shown. We conclude that the errors on the scaling violation data obtained 
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from the combination of MARK II and ALEPH data are too large to exclude any of 
our sets. 




X 

Figure 35: As Fig. but zoomed and without the data, for clarity. In addition to the NLO 
(solid) and LO (dotted) results for our central choice A^^^ = 350 MeV, the range covered by 
the NLO variable-A sets is also shown as dash-dotted curves. The upper one corresponds to 
A(^) = 150 MeV and the lower one to A^*^) = 400 MeV. 



5.3 The Gluon and Charm Content of the Photon 

Before the advent of HERA, information on the parton content of the photon came 
almost exclusively from the measurement of F2 in single tagged 77 collisions'^. To- 
gether with the assumption of vector meson dominance (VMD), the quark PDF's of the 
photon have been reasonably well determined from these data. The gluon, on the other 
hand, does not contribute to F2 at LO and is hence poorly known. Consequently, the 
quark distributions of the various parametrizations, e.g. in ||110| , |121| , |126|| , are similar. 



whereas the gluon distribution depends crucially on the choice of PDF's. 

As argued in subsection 5.1, the IPP process at the NLC in the LASER mode may 
be used to measure the gluon content of the photon more precisely. However, we will 



16 



In this work we are concerned only with the real photon. For a comparison of real and virtual 



photon PDF's, see [125| 
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have to wait several years for the first data from the NLC - if it will be operated in 
the LASER mode at all. 




Figure 36: The impact of the gluon in the photon on 77-spectra of charged particles in 
photoproduction. The NLO results [|116|| for ^ = 0.6 with our set |13] of charged hadron 



FF's are compared to preliminary data from HI The cross section has been integrated 

over pj- > 3 GeV. The solid line gives the result with the GRV photon PDF's [110|, whereas 
the dot-dashed line is the result obtained with the GS photon PDF's [121| . The dotted 
(dashed) line results from switching off the gluon in the fragmentation (in the PDF's). The 
dot-dash~dotted curve is discussed in the text. 



In the meantime, HERA is the ideal place to learn more about the gluon content 
of the photon. In particular the rapidity spectra show a significant sensitivity to the 
gluon PDF of the photon. Fig. ^ shows the preliminary HI data on inclusive charged 
hadron production in photoproduction, integrated over pr > 3 GeV, as a function of 
the rapidity ||113|| . The prediction of our NLO analysis with set |]13[ of charged hadron 
FF's is plotted as solid line. As in the comparison with the px spectrum in section 
5.1, we use the GRV fTTUl photon PDF's and the CTEQ3M proton PDF's and 
choose the optimized scale ^ = 0.6 in our calculation. The agreement with the data 
is near perfect over the whole measured region of — 1 < 2/iab < 1- When using the 
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GS ||121|| PDF's instead (dot-dashed) the theoretical prediction decreases only slightly 
and is still consistent with the data. However, both sets are constructed with similar 
prejudice concerning the VMD input | |1 25| ] and it would therefore be misleading to 



take this small difference as a measure of the sensitivity of the observable to the gluon 
density in the photon. The sensitivity is revealed much more clearly by switching off 
the gluon in the photon, which gives the dashed curve. It is clearly inconsistent with 
the data in both normalization and shape. Moreover, the rapidity spectrum is sensitive 
to the gluon density at both small and large x^, where x-y is the fraction of the photons 
momentum that is carried by the gluon. A toy distribution that is similar to that of 
GRV around = 0.1 but zero else results in the dot-dash-dotted curve. This curve 
approaches the GRV result in the forward region (positive y\ah) and is identical to the 
no gluon result in the backward region. 

We conclude that the gluon distribution can be probed by IPP in photoproduction 
from ~ 0.05 to x^ considerably above 0.2. However, the rapidity spectrum is even 
more sensitive to the gluon FF. The result of turning off the gluon fragmentation (in 
the cross section) is plotted as dotted curve. Obviously, the precision of any future 
measurement of the gluon content of the photon in this process depends heavily on the 
uncertainty of the gluon FF. This holds in particular in the forward region where the 
gluon fragmentation contributes up to 80% of the cross section. Our improved charged 
hadron set should be sufficiently precise to permit an analysis along these lines. 

The charm distribution in the photon is not well constrained by the F2 data, either. 
At HERA, it can be probed by inclusive production of charmed mesons in photopro- 
duction. Precise data have become available for the production of D*^ mesons, recently 
r7| , pU[ . The FF's are well known as has been demonstrated in section 5.1. The main 



limitation for the measurement is the experimental error. Fig. shows the compari- 



son of our NLO predictions with data from ZEUS pO] for the rapidity spectra of the 



cross section integrated over > 3 GeV. The solid curve in Fig. |37| a represents our 



standard choice of the GRV set of photon PDF's, the CTEQ4M |TT3] proton PDF's, 
and scale ^ = 1. For this choice, we find good agreement with the measurement in the 
central and backward region, whereas it is about I.80" low for the bin in the forward 
region. The dotted curve represents the direct contribution. It is far below the data. 
When only the charm distribution in the photon is switched off, the dashed curve re- 
sults. It almost coincides with the dotted curve, except for the forward region, where 
it is slightly higher. Clearly, the dominant resolved part of the cross section is almost 
entirely due to the charm density. Therefore, a measurement of the rapidity spectrum 
of D*^ corresponds, after subtraction of the well determined direct part, directly to 
the measurement of the charm density in the photon. 

Fig. P7| b shows the NLO predictions of various PDF sets, compared to the same 
data. In addition to GRV, the set of GS ||121|| is plotted as dotted line. ACFGP 
presents two sets, with a massless charm (dashed) and with a massive one (dash- 
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Figure 37: The impact of the charm quark in the photon on ry-spectra of D*^ production 
in photoproduction. Our NLO predictions are compared to data from ZEUS [20|. The sohd 
Une, calculated with the GRV photon PDF's [110| gives reasonable agreement with the data, 
a) The result with no charm in the photon and the direct contribution are also given, b) 
Sensitivity of the cross section to the differences between the photon PDF's. The NLO results 
are shown for the sets of GRV GS [|l2|], and ACFGP |26|. 



dotted). GS is slightly favored over GRV by the data, whereas the massive set of 
ACFGP compares least favorable to the data. Unfortunately, the statistics are not 
sufficient at the moment, to exclude any of the above sets. 
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6 Summary and Outlook 



In this thesis, we set out to extract fragmentation functions for charged hadrons to 
NLO accuracy from e^e~ annihilation data solely. Former fits revealed that it is not 
possible to unambiguously determine the gluon FF with data on unpolarized cross 
sections alone We thus used data on the longitudinal polarized cross 



section to constrain the gluon fragmentation. 

To this end, the NLO longitudinal coefficient functions are required. We completed 
the calculation of P.J. Rijken and W.L. van Neerven |]^ by the color class CpT^ and 
checked the pure singlet and the class NpTuCp of the nonsinglet. These coefficient 
functions were then used together with those for the transverse cross section, to fit 



FF's to ALEPH data on the longitudinal and the unpolarized cross section. In 
the unpolarized case, ALEPH distinguished three different data samples for (i) mainly 
M, ci, s quarks (ii) mainly c fragmentation and (iii) a h quark enriched sample. Our 
analysis yields well constrained light quark, c, and h quark fragmentation functions for 
charged hadrons. Owing to the inclusion of the NLO longitudinal coefficient functions 
in the fit, the gluon is also well constrained in the important small x region. However, 
it is not constrained well for intermediate and large x. As the gluon fragmentation is 
effectively probed at relatively small x in many processes, this analysis nevertheless 
represents a significant improvement over our earlier sets. 

In particular, the new sets could be used to measure via inclusive particle pro- 
duction in photoproduction, DIS or pp scattering, where the gluon dominates the cross 
section. For a consistent analysis along those lines, the value of A must be varied not 
only in the formula for the strong coupling, but also in the PDF's and the FF's. It is 
thus necessary to perform separate fits for a range of specified A values. For the PDF's 
of the proton this has been done some time ago by A. Vogt ||124|| . In this thesis, we 



provide an analogous set of functions for the fragmentation process, for the first time. 

Our fragmentation functions may also be put to good use in constraining the charm 
and the gluon densities of the resolved photon from rapidity spectra as measured at 
HERA. With the precision that is expected from the luminosity accumulated in the 
past years, this is arguably the best place to learn about the charm PDF. Finally, these 
FF's will give particularly precise predictions for cross sections at future experiments 
like LEP2 and NLC First data from LEPL5 |p.27|| are already awaiting comparison 



with theory. All of this has hardly been tackled, yet, promising fascinating insights to 
be gained in future work. 

Another application of fragmentation functions is to study scaling violation in e^e~ 
annihilation. In we fitted A, along with the other parameters, to data at 29 GeV 
and at 91.2 GeV. We obtained very sensible values for asijnz) of 0.118 (0.122) in 
NLO (LO), however, we could not assess their uncertainty by this simple approach. A 
dedicated study with complete error analysis can give competitive errors for the value 
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of the strong coupling, as demonstrated, e.g., by ALEPH |45| . 

Universality, a key prediction of the factorization theorem, can also be tested by 
applying FF's that were fitted to e~^e~ annihilation data to other processes. This has 
also been demonstrated successfully for ep-, pp-, and 77-scattering, and is briefly 
reviewed in this work. With the improved set of FF's, it will be possible to reduce the 
theoretical uncertainties considerably. 
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A The Strong Coupling and Aqcd 



The coupling constant of QCD satisfies the renormahzation group equation 



2dOis ni \ 



(A.l) 



This differential equation determines the scale dependence, hence is known up to a 
constant. To NLO, the first two terms in the expansion are kept, 



a: 



47r yvKy 



33 - 2Nf 



306 - 38A^F 



3 3 
and to this order, the solution of eq. ( |A.1| ) can be expanded to 

2 



(A.2) 



27r 



/3oln(/iVA2) 



A Inln(^VA^) 
i5l ln(/iVA2) 



(A.3) 



up to terms of order [ln(/i^/A^)]~'^ that can not be consistently included at NLO. 
(Other approximations are in use also, introducing numerically small ambiguities in 
the interpretation of as and A — see ||128|) To LO, the logarithmically suppressed 



term in the square brackets is absent. The parameter A is the one constant, mentioned 
above that fixes the boundary conditions of ( A.l| ). It can not be calculated from 
perturbative QCD but must be determined from experiment. 

In this work, A is used as a phenomenological parameter to fix the value of the 
strong coupling at some point. For ease of notation this parameter is at places given 
without indices. However it does depend on a number of details in the calculation of 
the cross section. Great care should be taken when comparing A values of different 
analyses. First, this QCD parameter depends on the renormahzation scheme used. 
In this work the MS-scheme pO[ is used throughout. For two distinct schemes, say 



scheme A and scheme 5, the coupling constant is identical to LO but differs by a finite 



renormahzation at NLO [ 128| 



of 
27r 



27r 



27r VV27r/ 



(A.4) 



The first two coefficients of the QCD /3-function are unchanged by such a transforma- 
tion so that /3o, /?! and C\ determine any scheme transformation (to NLO). In fact, only 
the leading coefficient of the /? function is needed for the transformation of A [ |128[j . 



A^ = A^ exp 



(A.5) 
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Second, as (3 depends on the number of flavors, it is clear that A does, too. When an 
effective theory is applied as done in this work, all the quark masses are taken to be 
massless in the calculation of the hard processes and the masses of the heavy quarks 
(c and b) are incorporated by the introduction of thresholds, only. This leads to a 
variable number of flavors, depending on the scale, and consequently the A parameter 
must be adjusted accordingly. The matching at the flavor-thresholds must be done in 
such a way that the strong coupling is continuous. The respective equations are rather 
intricate, for practical purposes it is sufficient to use the approximations presented by 
Marciano ||129|| , 
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and 



(A.6) 



where the parameters /i4 and /is specify the threshold scales for the transitions of 3 to 
4, respectively 4 to 5 flavors in the VFNS. In this work, the thresholds are set at twice 
the heavy quark masses: fi4 = 2mc ~ 3 GeV and /is = 2mf, ^ 10 GeV. Equations ( |A.6| ) 
apply to the next-to-leading A, in LO the factor in square brackets is absent and the 
relations are then exact. 
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Figure 38: The strong coupling with the Marciano scheme for A-matching [129]. Because 
of the limited accuracy of ( [A.6| ) , the ratio of the two-loop to the one-loop exhibits small 
discontinuities at the thresholds. 
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The accuracy of the A matching prescriptions ( |A.6| ) is demonstrated in Fig. |3^. The 
quantity plotted is the ratio of the two-loop (eq. (|A.3| )) over the one-loop as a 
function of the scale /x. As the matching prescription is exact in LO, the discontinuities 
at the thresholds are due to the error of the approximation for the NLO matching ( |A.6| ). 
The error is small compared to the effect of the higher order term in ( |A.3|) . 



p^{i),NLO 


^{4),LO 


p^{b),NLO 


^(5),LO 


as(A, mz) 


150 


52 


90 


33.5 


0.1036 


200 


71 


123 


46.5 


0.1081 


250 


90 


158 


60.0 


0.1119 


300 


109 


193 


74.0 


0.1152 


350 


129 


229 


89.0 


0.1182 


400 


149 


265 


104.0 


0.1209 



Table 18: Matching of A^'') to A^^). The A values are given in units of MeV. 



In photoproduction and in other processes with relatively low hard scales, only 
four flavors are active over much of the kinematical region and one commonly cites 
the A value for four flavors. In our new charged hadron fit we cite the A*^^-* values. 
They are chosen as in |p.24|| so that our set of FF's can be easily matched with the 
respective PDF's for a determination of a^. Table ^ lists the corresponding A^^^ 
values according to Marcianos matching prescription adopted here. The LO values are 
fixed by the requirement that the physical ^^(A^^^ m^) should not depend on the order 
of perturbation theory. 
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B Phase Space Integrals 



In order to perform the integration over four of the five degrees of freedom of the four- 
particle phase space, we regard the process in a specific system ||130|| . This is the CM 
system of partons 1 and 2, where the coordinate system has been rotated in such a 
way that of the five invariants in the irreducible set, ?/i2, yus, x^, 1/234, and 1/23, only 
2/23 depends on the angle cj). Also, only 2/234 and 1/23 depend on the second angle, 9. 
In this system, the angular integrations are relatively easy. With the definition of the 
integration measures 

l0[/(0)] = r^sin-2^0/(0) ^B-L) 
Jo 

= / -^v-%l-v)-^f{v) (B.2) 

JO I\v 

the phase space integrals required for the angular integrations of PS^^-* are given by 



U [1] = 1 



I.e. 



,2.„r(l-2e) 



Z TT 



COS 

I 

I 



2fc+l 







e INn 



COS 



COS 



(B.3) 
(B.4) 
(B.5) 
(B.6) 



for the integration over the polar angle 0. The normalization of eq. ( B.3 ) is determined 
from 



/ dzz^{l- zf = B{l + a,l + b) 
Jo 



r(i + a)r(i + &) 

r(2 + a + b) 



(B.7) 



The results for the azimuthal integration are (with the substitution v = 1(1 — cos 6)) 



I. [1] = 1 



i.e. A^^ = B(l -£,1 -e) 



B(l + k-eA-e) 



1 



v + D 
1 



B{l-e,l-e 
D + 1 



In 



D 



e 



Li2 ( ] - U2 



D 



D+1 



D(l + D) 



e{2D + 1) In [^^^] + 0{e''' 



(B.8) 
(B.9) 

0{e^) (B.IO) 

(B.ll) 
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where the expansion of Euler's Beta function in ( |B.8[ ) can be performed with the use 
of 



r(l + a) = ar(a) , r(i) = 1 

T{l + e) = l-^^e + Ull + a'2))e' + 0{e^) ; 7e~ 0.577215. 



The remaining azimuthal integrals are of the form 



with /c,/ e No 



(B.12) 
(B.13) 



(B.14) 



and are solved by recursive application of the decomposition 



(v + D) 



[v + oy 



D- 



v + D) 



(B.15) 



The above angular integrals are in agreement with | |131|| and have also been checked 
numerically. 

The integrals required for the energy integrations are too numerous to list them 
here. For the regular terms, the limit e — can be taken in the integration measures 
( p. loop and (|3.103|) and one has then to deal with standard definite integrals, merely. 
Those can be found in ||132|| . In the singular terms ( p.l06| ) the poles are extracted with 
the substitution 

z-'-' = --5iz) + (-) +Oie). (B.16) 



Due to the functional form of g{z,x,r) in ( p.l06| ), the plus distribution gives no con- 
tribution to the integral and up to 0{e'^) the singular term can be replaced by —1/e. 
However, the second energy integration must keep track of the terms proportional to 
e that, together with the pole, give rise to finite terms. For this purpose, the expres- 
sions in the second integral which is free of singularities are expanded in powers of e 
according to (|3.115| ). This, again, leads to standard integrals. 

Bookkeeping is facilitated by employing an algebraic computer program. In this 
work, MAPLE [|133|| has been used for the four integrations over phase space. 
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C Numerical Techniques for Evolution 



In this appendix, the different methods for the numerical solution of the Altarelli-Parisi 
evolution equations are discussed. Before turning to the details of the computer codes, 
the Mellin transform will be introduced as the basic tool for the more elegant of the 
two computer codes - the n-space code. 

The Mellin transform technique, however, is not just a convenient way to solve the 
AP evolution equations. The essential feature of the Mellin transformation 

rl 



/ 



n] = 



is that it renders convolutions to products. 



n] 



f 



n 



(C.l) 

(C.2) 
it will be 



This can be employed for many calculational purposes, see |p.34|| . Here 
exploited to reduce the system of differential equations describing the evolution of 
FF's to a system of linear algebraic equations. This system is solved analytically ||135 



and the solution is then transformed back to x-space via the inverse Mellin transform 

f{x) = ^J^dnx-^f'\n). (C.3) 

The Mellin technique has the additional advantage that one keeps exactly those terms 
that have to be kept in the order in perturbation theory considered, whereas in x-space 
evolution one can not help but sum up part of the higher orders. The x-space evolution 
is in that sense somewhat inconsistent. 

A drawback of the Mellin technique is, however that the input for the evolution 
must be known analytically in order to make use of the techniques advantages. Even 
if it is known analytically, it is not always possible to perform the Mellin transform 
analytically which then spoils the techniques main numerical advantages. For details 
of the numeric in n-space, see |p.36| , p.37| . 



We will closely follow | ]135| | in the ensuing sketch of the analytical solution in n- 
space. With the introduction of the new variable 



t 




(C.4) 



the set of Mellin-transformed AP equations (|2.21|) - ( 2.24 ) can be written as 



dtMc 



7(+)^(+),i 



.gMg 
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The convolutions have been rendered to products of the new functions M = Z^I^^l and 
7 = p[^^]. For the standard parametrization ( [4.1|) , the moments of the FF's are Beta- 
functions ( IB.?! ). The timehke anomalous dimensions 7 can be found in | ]138| | . The 
introduction of matrix notation 



7 



7s Iq^G 



M 



Mg 



(C.5) 



further simplifies above equations is disposable as it does not enter the cross 

sections (|3l5|)). 



aiM(+),i = 7(+)M(+), 
dtM = 7 M . 



(C.6) 
(C.7) 



In the following, matrix notation is implicit where appropriate. The solution of the 
evolution equations can be expressed by the action of evolution operators E 



M{t,n) =EM{0,n) 



where for the non-singlet 



E{+){t,n) = 
and for the singlet 
E{t,n) = 



1 + ^s{t)-as{0) ( A (0) . . _ 1 (1) 



ei + 



27r 



as{t) - as{0) 



^7(+)H-^7wH 



(C.9) 



+ 



62 + 



27r V 

asjt) - asjO) 
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2 „ , e2Rei 
-Jo'"'''' + X.-X.-hPo, 
( 2 e^Re2_ 



The Matrix R is defined as 



The Aj are the eigenvalues of the matrix 7*^°-' 
1 



Ai,: 



and the matrices 



G 



1 



ei,2 



Al,2 — ^2,1 
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7 



(C.IO) 



(c.ii: 



(C.12) 



(C.13) 



are projectors, with the properties 

6162 = , ei + 62 = 1 , 7^°^ = AiCi + A2e2 



p2 



ei,2 



(C.14) 



Here, 1 stands for the one in the matrix-space. The evolved FF's in x-space are then 
obtained via numerical back-transformation of the M{t,n), according to ( |C.3| ). 

In the spacelike number of computer codes for the numerical solution of the 

evolution equations have become available by now [|139 



The theoretical uncertainties 
and the differences between the various codes are discussed in ||140|| . In the timelike 
case, progress has been quite recent. Nason, e.g., wrote the possibly first evolution 



program for the timelike region [^. Spira and Kniehl completed computer codes in 
the recent past |141|, in addition to the codes by the author that have been employed 
in p!3|,p!^,p!5|,|3T|,|142|| . Our codes for both the n-space and the x-space will be briefly 



discussed below. 

(i) n-space: In n-space, the AP-equations are solved analytically, as outlined in this 
section. An adaptive routine is used for the back-transformation into x-space 



The contour is indicated in Fig. |3^ as the dotted line. It runs parallel to the imaginary 
axis, through (7,0), and is closed at infinity. The real part of the kernel in (|C.3|) is 
symmetric to the real axis so that it is sufficient to integrate over half the contour. For 
large imaginary parts of n, the kernel rapidly approaches zero so that the integration 
can be terminated at rather small values of |max| = (9(20). Both |max| and 7 are 
chosen adaptively. The precision can thus be set to any desired value, usually it is of 
the order of 0.5%. For the calculation of the complex logarithmic gamma functions the 
NAG FORTRAN library | ]143[ is used. The code has been checked numerically with the 
program by Spira ||144|| . 



Im 



max 



2^3 



Re 



Figure 39: Contour for the Mellin back-transformation, 
(ii) x-space: In this brute-force approach, the AP-equations are solved by iterative 



application of their integral form ( ^.341 ) where the convolution P®D is evaluated with 
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linear interpolation between n = (9(200) regularly spaced values and ( p.34[ ) is iterated 
of the order of 50 times. The accuracy of this routine is of the order of 1% for x < 0.8. 
We find perfect agreement with the program by Kniehl |JL44| . 



For both the x- and the n- space evolution, the accuracy deteriorates as x — > 1. 
As this is true also for the data (because a —>■ 0, too, and hence also the statistics), 
the impact of these numerical errors is small in most analyses. If desired, arbitrary 
precision can of course be obtained for any value of the momentum fraction x, albeit 
at high CPU time cost. 

A fundamental difference between the x-space and the ?7.-space techniques turns 
up in the way the perturbative expansion is truncated. Whereas one consistently 
keeps terms of C(a^) {0{as)) only in the NLO (LO) n-space evolution, the repeated 
evaluation of ( |2.34| ) in the x-space approach amounts to a partial resummation of 
higher order terms. However, as the deviation is numerically rather minor it may be 
disregarded for most practical purposes. Fig. ^ shows the ratio of the FF's obtained 
at mz when using the n-space technique over those obtained in the x-space. In both 
cases, our new NLO set for charged hadrons was used as input at /iq. The largest 
deviations from 1 show up in the fragmentation of the gluon (solid line). The light 
quarks (dashed) are also affected slightly whereas the charm and bottom FF's are 
insensitive to the choice of technique. In general, the effect is largest for the partons 
that dominate the evolution. (See [p.40|| for a detailed discussion of the ambiguities in 
the closely related spacelike case.) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 40: Ambiguities of the solution of the AP equation. The curves show the ratios of the 
results of the n-space evolution over those for the x-space evolution. Our new NLO sets of 
FF were evolved from fiQ to mz- The solid, dashed, dotted and dash-dotted lines correspond 
to the fragmentation functions of the gluon, the light quarks, the c and the b, respectively. 
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D Coefficient Functions 



The order a° coefficient functions in eqns. ( |3.25| ) - ( |3.28| ) are given by 



(D.l) 



<CtT'^'\z) = else. 



(D.2) 



Only at 0{as) do the coefficient functions differ from the naive parton model expecta- 
tion [Hi! 



^-{l-z)+ 4^ - (1 + ^) (ln(l -z) + 2 \n{z)) 



+ 2 



ln(l -z) 
1 - z 



z) = €lf(z) = Q 



1 + 



= Cf[1] 



(In(l-z) +21n(2)) -2 



I- z 



(D.3) 
(D.4) 

(D.5) 
(D.6) 



1 - z 



(D.7) 
(D.8) 



For the factorization, also the a functions proportional to e as introduced in ( |3.107| ), 
(|3.109|) are presented, 



ail{z)=CF[l-2\n{z)-\n{l-z)] 



A- (2-2 ln(z) - ln(l - z)) 

z 



(D.9) 
(D.IO) 



The C^"* are not required in this work, see P] for the expressions. The same applies 



to the and C^'^ that can be found in [pOl. All of the CV are collected in the 



.(2) 



.(2) 



following equations for convenience. The color class CpT^ of the longitudinal cross 
section, specified by the coefficient function C£ is presented here for the ffist time. 
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Class NfTrCf of C^^ ^ as well as C^ ^ have been checked and the result of [0], where 
all the other classes are taken from, has been confirmed. 



= C|{45i,2(l-z)+85i,2(-2)-12Li3(-z)+81n(l + z)Li2(-z) 
+41n(z)Li2(-z) + 4C(2) ln(l + z) - 4C(2) ln(l - z) 
+41n(z) ln2(l + z)-2 \n^{z) ln(l + z) - 4C(3) 

+ (^-4 + y ^ - 8z - ^z^^ (Li2(-z) + Hz) ln(l + z)) 
-3Li2(l -z)+ \n{z) ln(l - z) + (^4 - 8z - C(2) 

+ ln2(l - ^) + (^-^ + 4z + ^^3^ H{z) + + ^) - ^) 

/17 12 1 6 8 9\ , , , 147 12 1 9 8 , 

+ V -z+ -z^] \niz) \ z + -z^ 

VlO b z b b ) ^ ' 20 5zl0 5 

+A^cCf {-25i,2(l -z)- 45i,2(-z) + 6Li3(-z) 

-41n(l + z)U2{-z) - 2C(2) ln(l + z) + 2C(2) ln(l - z) 

-21n(2)Li2(-z) - 2hi(z) ln2(l + z) + \n^{z) ln(l + z) + 2C(3) 

6 1 4 

5;^ + ^^+5^ 



+ ( 2 - -— + 4z + -z^ ) (Li2(-z) + ln(z) ln(l + z)) 



+C(2) (^4z + - {iz + \n\z) - ^ ln(l - z) 

f Ti 6 1 2 4 2\ , . ^ 1729 6 1 49 4 

+ ^ h -z z^ Infz) H z z 

V 30 5z 5 b J ^ ' 180 5z 30 5 

r 2 25 2 1 

+iV^T^C^ |-(ln(l - z) +ln(z)) " y + 3^ j (D-H) 

r OS 1 f^o 

>(z) = A^i^r^^C^ 4Li2(l-z)+41n(z)ln(l-z)+61n2(z)- — -4- + — z-4z2 

+ (^^ - 4z + ^z^) ln(l - ^) + (^-8 + - 8z + ^z^) ln(z)| (D.12) 



.^.{2)^.^ _ r^„^„J^T^.n _^)_^(2)) (_4 + i6i-i64) - -z^ 

\ Z Z J 

'2 4 1 \ ^ , . /8 32 o 16 1 



z z^/ 5 
+Li2(z) ( ^ + 7^ ) + Li2(-z) f ^ - ^z^ + 



3 

3 ' 5zV ' "^^^ V3 IS'^ ' 5 z2 
, / 1 1 \ , , , / 16 22 44 1 

Sw-, N 32 2 32 , , 16 1 ,^ , 16 , ,\ 
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641 56 8 32 ol ^ 

+ z + — zH (D.13) 

5 z 5 15 15 J ^ ^ 
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C%{z) = C| { (-^ + ^1 + ^^3^ (Li2(-z) + Hz) HI + z)) 

+4Li2(l - z) + 41n(z) ln(l - z) + -|c(2)2^ 

15 

+ ^z^^ H{z) + (^-6 + 8^ - ln(l - z) 

( 2 24 1 56 8 , . N 6 24 1 62 80 

+ \ z z^ HA H \ -2 

V 5 5 z 15 15 / ^ ^ 5 5 z 15 15 

^NcCf I (^8 + 8^^ (Li2(-2) + HA + 2)) - 16^Li2(l - z) 

-161n(z)ln(l - z) + C(2) f-16 + 24-) + (^-4 + 4^) H{\-z) 



1\ , / 116 1 4 ,\ , 

12 + 16- Hiz) + 36 - — — + 4z - -z^ In 1 
2; / V 32; 3 / 

1 . „ 4^2V ' ^ 80 . 112 1 40 



+ , 28 - -- + 8. - ln(.) -_ + _--_. + -.^ J (D.14) 



The plus-distributions that appear in the coefficient functions and in the sphtting 
functions can be defined via the convolution with a regular, so-called test-, function. 

1 1 
^dz \]{z)\^ g{z) ^ jdz f{z) {g{z) - g{l)) (D.15) 


The subtraction is to be taken at the locus of the singularity, which in this work is 1, 
mostly. The convolution of a plus-distribution with a regular function is then given by 
adding the missing contribution to eq. ( p.l5| ). 



g{x) / dzf{z) . 
Jo 



(D.16) 



Nielsen's generalized polylogarithms ||146|| are defined as 

Li2(z) 



1 

113(2;) = f dy 
Jo 



dy- 
y 

1 p ^ ln2(l - y) 



Li2(y) 



Riemann's zeta function has the values 



TT 



((2) = — ^ 1.644 934 
6 

C(3) ^ 1.202 057, 



(D.17) 
(D.18) 
(D.19) 

(D.20) 
(D.21) 
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